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SECTION  1 


INTRODUCTION 


1-1.  BACKGROUND 

The  MISERS  BLUFF  II  (MB  1 1 )  Cloud  Sampling  Program  was  initiated 
to  provide  the  empirical  data  required  to  determine  how  dust  cloud 
characteristics  of  single  burst  environments  interact  to  form 
multiple  burst  environments.  State-of-the-art  dust  cloud  modeling 
relies  on  the  principle  of  linear  superposition  in  defining  the 
flyout  threats  for  advanced  basing  concepts,  hence,  validation  of 
this  hypothesis  has  become  a  very  important  consideration. 

The  high  explosive  (HE)  tests  in  the  MB  1 1  series  (FY78)  provided 
a  unique  opportunity  to  gather  the  dust  particle  environmental  data 
from  both  a  single  burst  dust  cloud  and  a  dust  cloud  generated  by  six, 
simultaneous,  interacting  bursts .  Science  Appl ications ,  Inc.  (SAI), 
under  the  sponsorship  of  the  Aerospace  Systems  Division  of  the  Defense 
Nuclear  Agency  (DNA/SPAS),  directed  and  performed  extensive  dust 
particle  sampling  experiments  during  both  events  in  the  MB1I  series. 

The  experiment  rationale  and  procedures  utilized  in  sampling  the  dust 
clouds,  along  with  the  preliminary  results,  have  been  documented  by 
SAI  in  DNA  4729F,  "MISERS  BLUFF  II  Cloud  Sampling  Program:  Procedures 
and  Preliminary  Results"  (Ref.  1).  The  purpose  of  this  document  is 
to  summarize  the  extensive  analysis  that  has  been  conducted  with  the 
empirical  data  gathered  during  the  experimental  phase  of  the  program. 

1-2.  DATA  BASE  FOR  THE  MB  1 1  DUST  CLOUD  ANALYSIS 

The  majority  of  the  dust  cloud  analyses  focused  on  the  in-situ 
dust  particle  spectral  data  obtained  from  both  dust  clouds  by  Particle 
Measuring  Systems,  Inc.  (PMS).  This  data  was  gathered  from  an  airborne 
platform  (Cessna  206)  as  it  penetrated  the  single  and  multiple  burst 
dust  clouds  with  approximately  20  sampling  passes  through  each  cloud. 

The  PMS  data  consists  of  dust  particle  number  concentrations  for  the  size 


range  from  %3-750  pm  (single  burst  cloud)  and  from  '0.1-3000  inn  (mul¬ 
tiple  burst  cloud)  with  M20  meter  (400  foot)  horizontal  resolution 
within  the  clouds.  The  results  of  this  effort  and  a  description  of 
the  particle  sizing  spectrometers  have  been  documented  by  PMS  in 
DNA  4951F,  "Results  of  the  MISERS  BLUFF  II  Aircraft  Dust  Particle 
Sampling  Experiments"  (Ref.  2). 

Supplemental  data  sources  utilized  in  this  analysis  includes 
data  obtained  by:  (1)  Science  Applications,  Inc.  (SAI)  -  large  particle 
fallout  from  the  multiple  burst  cloud  at  times  prior  to,  and 
coincident  with,  airborne  sampling;  (2)  Air  Force  Technical  Applica¬ 
tion  Center  (AFTAC)  -  filter  paper  samples  containing  dust  from  the 
top  of  the  MBI I  multiple  burst  cloud;*  (3)  Technology  International 
Corporation  (TIC)  -  high  resolution  cloud  photoqraphv  for  thirty 
minutes  of  both  events  in  the  MBII  series;  and  (4)  an  Air  Force  6th 
Weather  Squadron  (Mobil)  Unit  -  detailed  meteorological  data  from  the 
surface  to  cloud  top  (approximately  4600  meters  ( 1 5kf t )  above  ground 
level  (AGL))  for  each  event  in  the  series. 

In  all,  it  is  felt  that  the  data  base  for  the  MISERS  BLUFF  II 
dust  cloud  analysis  is  the  most  extensive  and  complete  that  has  been 
obtained  to  date  during  any  dust  cloud  sampling  program. 

1-3  REPORT  FORMAT 

Due  to  the  immense  quantity  of  dust  particle  spectral  data 
obtained  from  both  MBII  dust  clouds,  a  major  effort  in  producing 
this  document  was  directed  at  formating  the  data  in  a  manner  such 
that  little  or  no  explanation  would  be  required  in  interpreting 
the  way  the  data  has  been  presented  or  conclusions  that  have  been 
drawn  from  it.  Data  sets  from  the  single  burst  and  multiple  burst 
dust  clouds  have  been  made  consistent  with  each  other,  i.e.,  the  same 


*AFTAC  also  fielded  aerosol  spectrometers  which  served  as  back-up 
instruments  to  the  LAS-X  instrument  fielded  by  PMS  for  the  multiple 
burst  event.  Since  the  LAS-X  instrument  functioned  normally,  the 
AFTAC  spectral  data  have  not  been  reduced  or  included  in  this  analysis. 


dust  particle  data  has  been  presented  for  both  clouds,  in  order  to 
enable  readers  to  make  their  own  comparisons  of  the  salient  features 
of  the  single  and  multiple  burst  dust  clouds.  Further,  the  desire  to 
have  a  "complete1'  document  has  led  to  the  inclusion  of  several  compre¬ 
hensive  appendices  to  supplement  the  material  presented  in  the  main 
body  of  the  report. 

Section  2  provides  a  detailed  summary  of  the  spectral  data 
obtained  from  the  single  burst  dust  cloud  (MB  1 1  - 1 )  and  also,  describes 
the  gross  mass  loading  features  of  the  cloud  at  ten  and  twenty  minutes 
after  detonation.  Section  3  provides  an  identical  description  of  the 
multiple  burst  dust  cloud  (MBII-2).  Section  4  is  devoted  to  detailed 
comparisons  between  the  mass  loading  features  of  each  cloud.  Section  5 
is  included  to  discuss  the  significant  uncertainties  in  calculating 
total  mass  lofted  by  each  detonation  and  how  these  uncertainties 
would  effect  the  calculations.  Section  6  contains  the  conclusions 
and  recommendations.  Appendices  A-D  are  included  to  provide  the 
reader  with  detailed  data  on  several  other  dust  cloud  parameters 
obtained  from  the  measured  spectral  data  and,  also,  to  provide 
"back-up"  material  for  Sections  2,  3,  4,  and  5. 


SECTION  2 


THE  MISERS  BLUFF  II  SINGLE  BURST  DUST  CLOUD 

2-1  MB I I - 1  IN-SITU  SAMPLING  DATA  SUMMARY 

Although  more  than  20  sampling  passes  were  flown  at  various 
altitudes  through  the  MB  1 1 - 1  dust  cloud,  only  those  passes  that  were 
made  while  the  aircraft  was  ascending  to  the  top  of  the  cloud  are  in¬ 
cluded  in  the  mass  loading  calculations.  Further,  sampling  passes  1 
and  2  were  accomplished  below  the  cloud  cap  and  outside  the  visible  stem, 
hence,  data  obtained  from  them  have  not  been  included  in  the  calculations. 
Similarly,  passes  10  and  11  were  performed  while  the  aircraft  was  maneu¬ 
vering  in  and  out  of  the  cloud  and  they  have  also  been  excluded. 

Table  2.1  summarizes  and  profiles  the  in-situ  sampling  mission 
accomplished  within  the  MB  1 1  - 1  dust  cloud.  Pass  numbers,  pass  mid¬ 
point  time-after-burst  (TAB),  and  sampling  altitude  are  tabulated  at 
the  top  of  the  table.  Pass  numbers  in  parentheses  and  the  altitudes 
of  the  top  and  bottom  of  the  cloud  layer  they  represent  are  enumerated 
on  the  left  side  of  the  table.  Pass-averaot d  mass  concentrations, 
as  calculated  from  the  measured  particle  number  concentration  (see 
Appendix  B  for  details),  are  given  for  dust  particle  sizes  -  47  .  m  and 
■47  ,m  for  each  pass.  The  four  columns  on  the  right  hand  side  of 
Table  2.1  contain  calculated  values  of  total  pass- averaged  mass  con¬ 
centration  (all  particle  sizes),  cloud  volume  represented  by  sample,  mass 
in  cloud  layer  sampled,  and  cumulative  mass  in  cloud  as  measured  from 
the  cloud  top.  Volume  calculations  have  been  made  assuming  that  the 
dust  cloud  is  circular  at  each  samplinq  level  and  that  the  length  of 
the  sampling  pass  represents  the  diameter  of  this  circle  (more  will 
be  said  of  this  assumption  in  Section  5).  All  number  entries  ir,  the 
Table  are  given  in  an  abbreviated  format,  e.g.,  7.94+8  -  7.94x10  . 

Several  trends  in  the  data  become  apparent  in  Table  2.1.  Tirst, 
the  mass  concentrations  are  generally  decreasinn  with  time  (increasing 
altitude  and  pass  number).  This  is  due  mainly  to  cloud  diffusion  under 
the  influence  of  the  significant  wind  shear  aloft  measured  at  detonation 
time.  This  wind  induced  cloud  diffusion  is  also  evidenced  n  t  tie  T 1 f 
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TABLE  2.1  MISERS  BLUFF  II-l  In-Situ  Sanpling  Data  Summary 
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photographic  coverage  of  the  MB  1 1 - 1  dust  cloud  (Ref.  5).  This  fact 
is  further  supported  by  noting  that  although  mass  concentrations  are 
decreasing  with  time,  there  is  a  correspondi ng  increase  in  total  mass 
at  altitude;  the  maximum  being  found  at  1280  m  (Pass  12). 

Another  very  interesting  and  important  feature  revealed  by  the 
data  is  that  the  mass  contribution  from  dust  particles  2.47  pm  comprises 
only  v,i5  percent  of  the  total  mass  lofted.  This  fact  was  unexpected 
based  on  a  preliminary  analysis  by  PMS  of  subsurface  soil  samples  taken 
from  ground  zero  and  sized  by  one  of  their  probes  (Ref.  2).  That 
analysis  indicated  a  mass  peak  at  approximately  200-300  pm.  It  is 
apparent  that  either  the  flow  fields  generated  by  the  MB  1 1 -1  fireball 
rise  were  not  sufficient  to  keep  many  particles  of  this  size  aloft 
for  sustained  periods  of  time  or  that  these  larger  particles  were 
actually  "silty"  aggregates  and  were  broken  up  by  the  high  "pressures" 
from  the  cratering  and  particle  ejection. 

Based  on  the  in-situ  sampling  data  obtained  from  the  MB  1 1  - 1 
dust  cloud  and  assumptions  concerning  cloud  geometry  and  soil  density, 

g 

the  total  mass  lofted  in  the  single  burst  event  was  7.9x10  grams 

(highlighted  value  in  Table  2.1).  This  number  compares  very  favorably 

8 

to  that  value  calculated  by  PMS  (8.0x10  g)  in  an  independent  analysis 
performed  on  the  MBII-1  dust  cloud.  (Preliminary  analysis  of  MB  1 1  - 1 
crater  and  ejecta  measurements  by  AFWL  and  CERF  have  shown  that  approxi- 

O 

mately  10x10  g  of  soil  were  unaccounted  for  after  the  burst.) 

2-2  MBII-1  DETAILED  SAMPLING  PASS  DATA 

During  each  sampling  pass  through  the  MISERS  BLUFF  1 1  - 1  dust 
cloud,  particle  number  concentrations  were  measured  with  two-second 
resolution  (400  foot)  within  the  cloud.  From  the  measured  number  con¬ 
centrations,  it  was  possible  to  calculate  mass  concentrations  every  two 
seconds  and  plot  time  histories  for  each  sampling  pass  tor.  idet ed  in 
the  analysis.  Figures  2.1  through  2.13  contain  the  detailed  sampling 
pass  data  derived  from  these  measurements .  With  the  exception  of 
sampling  passes  1  and  2  (data  not  included  in  the  analysis),  pass  data 
arepresented  in  pairs  of  figures;  rigure  2.X  (a)  contains  pass  averaged 
perticle  size  distributions  (PSD)  and  cumulative  mass  plots  while 
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Figure  2.X(b)  depicts  the  time  ("radial")  dependence  of  mass  concen¬ 
tration  within  the  dust  cloud. 

In  Figures  2.X(a),  the  particle  size  distributions  of  particles 
47  pin  are  plotted  with  circles  while  the  size  distributions  of  par¬ 
ticles  47  nm  are  plotted  with  squares.  A  straight  line  representing 
the  "best  fit"  has  been  drawn  through  the  large  particle  data  so  that 
the  slope  of  the  PSD  can  be  easily  determined.  The  dashed  line  at  the 
bottom  of  each  figure  is  merely  the  upper  decade  of  the  large  particle 
PSD  plotted  on  the  same  scale  as  the  small  particle  PSD.  The  cumula¬ 
tive  mass  versus  particle  diameter  plot  is  the  curve  that  runs  from 
the  lower  left  to  the  upper  right  in  Figures  2.X(a);  values  of  cumula¬ 
tive  mass  (percent)  are  given  on  the  inside  right  hand  scale. 

Figures  2 . X ( b )  depict  time  histories  of  mass  concentration 
within  the  MB  1 1 - 1  dust  cloud.  Each  figure  contains  the  contribution 
to  mass  concentration  from  particles  47  urn  (circles  connected  with 
solid  line)  and  particles  47  pm  (squares  connected  with  dashed  line). 
Values  of  total  mass  concentration  can  be  obtained  from  Appendix  B 
(Table  B.l).  The  values  of  mass  concentration  in  Figures  2.X(b)  are 
logrithmic  values  and  were  plotted  in  this  manner  for  scale  compression. 

_4  r  3 

If  values  of  mass  concentration  were  calculated  to  be  10  cj/m  ,  they 

-  00 

have  been  plotted  as  10  ,  i.e.,  zero.  Pass-averaged  values  of  mass 

concentration  (e  in  g/m  )  are  represented  by  the  horizontal  lines  drawn 
through  the  data  (  47  ,.m  -  solid  line;  47  pm  -  dashed  line).  To  obtain 
pass  time  after  detonation,  subtract  13:05:00  from  the  pass  start  time. 

t  .  3  GENERAL  COMMENTS  CONCERNING  THE  MB1I-1  DUST  CLOUD 

Every  attempt  has  been  made  to  present  the  data  in  Fiqures  2.1 
fhrouoh  2.17  as  unambiguous  as  possible,  hence,  only  general  comments 
concerning  tire  sampling  pass  data  will  be  made  in  this  section  of  the 
report . 

2-3.1  Particle  Size  Distributions 

From  the  plotted  particle  size  distributions  depicted  in 
F inures  ?.X(a)  it  appears  that  ,  tor  particles  47  .  in,  the  pass-averaged 
size  distributions  encountered  in  the  MB  I  I  - 1  dust  cloud  can  be  approxi¬ 
mated  by, 


1/ 


d  >.  47  pm 


^  (#/cm3*pm)=  kd‘q  , 

where  d  is  the  particle  diameter  in  microns  and 

3.5  £  q  £  4.0. 

It  should  be  reiterated  that  this  range  in  q  is  for  pass- averaged 
data  -  it  would  probably  vary  significantly  within  each  two-second 
data  acquisition  interval  in  each  pass.  Further,  poor  counting 
statistics  for  particles  £  200  pm  lead  to  significant  uncertainty  in 
q  for  the  tail  of  the  distributions. 

The  size  distribution  of  particles  <  47  pm,  however,  appears 
to  be  approximately  exponential  in  nature,  i.e., 

(#/cm3*pm)  =  ke”q<^,  d  <  47  pm 

The  spectral  data  obtained  from  the  multiple  burst  dust  cloud  for 
particle  sizes  0.1-6  pm  (see  Appendix  C)  indicates  a  general  increase 
in  number  concentration  with  decreasing  particle  diameter,  hence,  the 
possibility  of  the  aerosol  dust  particles  being  distributed  log-normally 
can  be  ruled  out.  This  is  mentioned  because  most  current  fallout  models 
assume  the  radioactive  dust  particles  are  log  normally  distributed. 

2-3.2  Cumulative  Mass  Distributions 

With  the  exception  of  Pass  3,  where  only  about  3  percent  of 
the  mass  can  be  attributed  to  particles  47  pm,  all  subsequent 

sampling  passes  reveal  that  approximately  85  percent  of  the  dust 
cloud's  mass  is  attributed  to  particles  47 .  Further  shown  in  Figures 
2.x(a),  is  that  the  mass  median  diameter  a-pears  to  lie  between  25  urn 
and  30  .  m. 

2-3.3  Mass  Concentration  Time  Histories 

As  might  be  expected,  the  highest  mass  concentrations 
encountered  in  the  MB  1 1  - 1  dust  cloud  were  in  the  stem  at  4-5  minutes 
after  burst  (Passes  3  and  4).  Average  concentrations  in  the  stem  were 
approximately  0.3  g/m  with  the  maximum  mass  concentration  being 
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1.3  g/m  encountered  during  pass  3  (first  pass  in  stem).  For  all 
passes  made  within  the  cloud  cap  (passes  6-15),  maximum  values  of  mass 
concentration  are  approximately  3-4  times  the  average  concentration 
(see  Appendix  B,  Table  B.l  for  maximum  values). 

Passes  1  and  2  were  made  below  the  cloud  cap  and  to  outside 
of  the  stem,  therefore,  the  mass  concentration  time  histories  from 
these  passes  consist  of  dust  particles  falling  out  of  the  cloud  cap. 
During  both  passes  there  appears  to  be  an  area  near  the  cloud  center 
that  is  nearly  void  of  "fallout".  This  observation  is  also  revealed 
in  passes  4  (stem),  5  (stem),  7  (cloud  cap),  and  9  (cloud  cap)  where 
there  are  two  areas  of  high  mass  concentration  separated  by  a  central 
region  of  low  mass  concentration.  There  is  still  some  evidence  of  this 
feature  in  passes  14  and  15,  however,  at  these  late  times  this  might 
be  due  to  wind  shear  "pulling"  the  cloud  into  several  pieces.  The 
fact  that  this  "hole"  in  the  cloud  is  not  evidenced  in  the  mass  con¬ 
centration  times  histories  from  passes  6  and  8  is  probably  due  to  the 
aircraft  not  passing  through  the  central  region  of  the  cloud.  (A 
comparison  of  the  sampling  pass  lengths  of  passes  6  and  8  with  pass 
lengths  of  passes  5,  7  and  9  reveal  that  the  maximum  dimension  of  the 
dust  cloud  was  not  sampled  during  passes  6  and  8.)  Hydrocode  (HULL)  flow 
field  calculations  that  were  performed  in  support  of  the  DICE  THROW  HE 
detonation  (600  ton  ANFO  -  hemispherically  capped  cylinder  charge  config¬ 
uration)  predict  a  reverse  flow  or  downdraft  in  the  center  of  the 
cloud  at  5  seconds  after  detonation  that  persists  for  approximately 
10  seconds.  This  downdraft  could  prevent  soil  from  initially  flowing 
into  the  central  core  region  of  the  cloud.  (Overhead  photography  of 
the  DICE  THROW  cloud  (Ref.  4)  did  reveal  this  mass  loading  phenomena 
at  T+15  seconds.)  Further,  the  lack  of  any  evidence  indicating  any 
toroidal  flow  in  the  MB  1 1  - 1  dust  cloud  as  it  rose,  would  suggest,  that 
any  flow  field(s)  surrounding  the  early  time  toroidal  ring  of  dust  were 
probably  too  weak  to  replenish  this  central  core  region  with  dust  as 
the  cloud  was  rising.  This  would  account  for  the  low  mass  concent  rat  ion 
in  this  region  at  later  times. 
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MISERS  BLUFF  II-l  PASS#1 
START  TIME  (MST):  13:07:! 
ALTITUDE  AGL  (km):  0-34 
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MISERS  BLUFF  II-l  PASS  #3 
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MISERS  BLUFF  II-l,  Pass  3,  Mass  Concentrati 
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Figure  2.4(a)  MISERS  BLUFF  1 1  - 1 ,  Pass  4,  Particle  Size  and 
Cumulative  Mass  Distributions 
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2.4(b)  MISERS  BLUFF  1 1 - 1 ,  Pass  4,  Mass  Concentration  Time  History 
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2.5(b)  MISERS  BLUFF  1 1 - 1 ,  Pass  5,  Mass  Concentration  Time  History 


Number  Concentration  (Particles/cm 


Number  Concentration  (Particles/cm' 


Number  Concentration  (Particles/cm 


1 


1  10  100  1000 


1 

t 


E 

i 


ro 

E 


QJ 

-Q 


Particle  Diameter  (um) 

Figure  2.7(a)  MISERS  BLUFF  II-l,  Pass  7,  Particle  Size  and 
Cumulative  Mass  Distributions 
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re  2.8(b)  MISERS  BLUFF  1 1 - 1 ,  Pass  8,  Mass  Concentration  Time  Hi 
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MISERS  BLUFF  II-l,  Pass  13,  Mass  Concentration  Time  History 
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Fiqure  2 . 1? ( a )  MISERS  BLUFF  1 1  - 1 ,  Pass  14,  Particle  Size  and 
Cumulative  Mass  Distributions 
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2-4  RECONSTRUCTING  THE  MISERS  BLUFF  II-l  DUST  CLOUD 

In  order  to  gain  more  insight  into  the  gross  mass  loading 
features  of  the  MBII  dust  clouds,  it  was  necessary  to  reconstruct 
the  clouds  at  various  times  after  burst.  Although  the  measured  mass 
concentrations  are  very  useful  in  identifying  radial  variations  in 
dust  densities  at  one  altitude  and  one  time,  they  do  not  give  the 
reader  "snapshots"  of  both  clouds  at  any  given  time  for  comparison  purposes. 

Based  on  cloud  dimensions  at  a  specified  time  after  burst, 
particle  fall  velocities,  and  in-situ  sampling  data,  it  is  possible 
to  approximate  the  mass  concentrations  within  various  layers  of  the 
cloud  at  any  specified  time  after  burst.  The  times  selected  for  the 
MBII  dust  cloud  analysis  were  ten  (T+10)  and  twenty  (T+20)  minutes 
after  detonation.  At  times  prior  to  this  (pre-stabilization)  it  was 
felt  that  flow  fields  within  the  dust  clouds  (cloud  buoyancy  and 
ambient  lapse  rate  operate  interacti vely)  would  be  so  uncertain  that 
any  analysis  based  on  them  would  not  be  meaningful. 

The  basic  methodology  followed  in  reconstructing  the  dust 
clouds  was  to  take  pass  totalled  mass  spectral  data  and,  depending 
on  what  time  the  data  was  acquired  with  reference  to  the  specified 
cloud  reconstruction  time  (either  before  or  after),  raise  or  lower 
the  particles  to  the  cloud  layer  in  which  they  would  be  at  the 
specified  time.  For  instance,  if  the  cloud  was  to  be  reconstructed 
at  T+10  minutes,  some  particles  that  were  sampled  at  T+5  minutes  may 
be  on  the  ground  at  T+10  minutes.  Similarly,  particles  sampled  at 
T+15  minutes  would  have  to  have  been  in  cloud  layers  above  the 
sampling  altitude  at  T+10  minutes. 

The  distance  that  the  various  sized  particles  were  raised  or 
lowered  depended  only  on  when  they  were  sampled  and  their  terminal 
fall  velocities.  Pass  mid-point,  times  were  used  to  specify  the  times 
when  the  pass  totalled  mass  spectral  data  were  acquired.  A  further 
simplification  to  the  calculations  was  made  possible  by  assuming 
that  dust  particles  47  ,m  sampled  in  a  cloud  layer  would  remain  in 
that  same  cloud  layer-  over  the  periods  rtf  time  considered  in 


reconstructing  the  clouds.  This  assumption  is  acceptable  when 
one  realizes  the  terminal  fall  velocity  of  a  47  ,.m  particle  is  only 
7  m/min  and,  at  most,  twenty  minutes  of  dust  particle  vertical  motion 
has  been  considered  in  the  analysis.  (See  Appendix  D  for  assumptions 
concerning  particle  terminal  fall  velocities  used  in  this  analysis). 

The  final  variable,  dust  cloud  dimensional  data  at  various  times-a f'ter- 
burst,  was  obtained  from  either  photographic  documentation  of  the 
MB  1 1  dust  clouds  (Ref.  5)  or  from  the  sampling  pass  lengths  at  times 
near  cloud  reconstruction  times. 

2-4.1  The  MB  1 1 - 1  Dust  Cloud  at  Ten  Minutes  After  Detonation 

Based  on  the  methodology  and  assumptions  outlined  in  the 
previous  paragraphs,  the  MBII  single  burst  dust  cloud  was  reconstructed 
at  ten  minutes  after  detonation;  the  resulting  mass  loading  (g/nf*) 
calculations  as  a  function  of  altitude  are  presented  in  Table  2.2.  The 
format  for  Table  2.2  is  almost  identical  to  that  used  in  Table  2.1. 

Pass  numbers,  times-after-burst,  and  sampling  altitudes  are  presented 
on  the  top  of  the  table,  while  the  top  and  bottom  of  cloud  layers 
sampled  in  a  particular  pass  are  located  on  the  left-hand  side  of  the 
table.  Tabulated  values  in  Table  2.2  for  particles  47  in  and  47  .  m 
at  the  various  cloud  layers  are  now  mass  values  (g)  as  opposed  to 
values  of  mass  concentrations  (g/iir)  in  Table  2.1.  The  four  columns 
on  the  right-hand  side  of  Table  2.2  yield  values  of  total  mass  at 
altitude  (obtained  by  summing  mass  contributions  from  all  particle 
sizes  at  a  given  altitude),  cumulative  mass  as  measured  from  the  cloud 
top,  volume  in  cloud  layer  at  TUO  minutes  (calculated  based  on  the 
estimated  T ♦ 1 D  minute  cloud  diameter  of  1700  m),  and  calculated 
values  of  average  mass  concentrations  at  altitude. 

Several  features  of  the  tabulated  data  should  be  pointed  out. 

As  previously  mentioned,  t  tie  mass  contribution  from  particles  with 
diameters  47  m  was  not  moved  out  of  the  cloud  layer  where  the  particles 
were  sampled;  hence,  it  is  not  "spread  out"  by  time  into  othei  flood 
layers,  further,  the  mass  spectral  data  obtained  dining  sampling 
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UFF  1 1 - 1  Dust  Cloud  Reconstructed  at  Ten  Minutes  After  Detonatioi 
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pass  8  from  particles  '47  ,.m  is  not  "spread-out"  into  other  cloud 
layers  due  to  the  proximity  of  the  data  acquisition  time  and  cloud 
reconstruction  time,  for  mass  spectral  data  obtained  in  subsequent 
sampling  passes,  the  mass  is  distributed  in  and  above  the  sampling 
altitude  while,  for  sampling  passes  made  prior  to  pass  8,  the  mass  is 
distributed  in  and  below  the  samplinq  altitude.  Finally,  pass  data 
from  passes  3,  4,  a no  5  have  been  averaged  in  calculating  the  average 
stem  mass  concentration. 

One  apparent  con  t  rad  i  c  t  ion  does  exist  in  t  he  data  in  Table  2.2 
and  that  is,  starting  with  pass  12,  the  distribution  of  mass  extends 
above  the  visible,  MB  1 1 - 1  cloud  top  (  2300m).  There  ace  two  possible 
explanations  for  this  inconsistency  in  the  data:  (1)  the  large*  dust 
particles  are  probably  not  spherical  in  shape  and,  hence,  their  terminal 
fall  velocities  have  been  overestimated  by  the  assumptions  made  in  this 
analysis  (see  Appendix  D);  and  (2)  the  larger  particles  may  be  organics 
(low  specific  gravity)  whose  mass  contribution  could  be  ignored. 
Whatever  the  reason,  the  mass  lofted  to  altitudes  greater  than  tire 
visible  cloud  top  comprises  only  3.6  percent  of  the  total  (loud  mass; 
therefore,  it  is  not  felt  that  this  inconsistency  creates  a  signitirant 
perturbation  in  the  analysis. 

2-4.2  The  MB 1 1 - 1  Dust  Cloud  at  Twenty  Minutes  After'  Detonation 

following  an  analysis  identical  to  that  used  in  recons t rue t  i  i  u 
the  MB  1 1  - 1  cloud  at  T+10  minutes,  tire  dust  cloud  was  also  ret  ons  t  mk  '  ed 
at  twenty  minutes  after  detonation.  Table  2.3  presents  the  mass  loading 
calculations  for  tire  T  +  2D  minute  cloud  and  is  identical  in  fnimat  to 
Table  2.2;  the  observed  cloud  diameter  is  now  44001,!. 

Due  t.o  the  relatively  late  time  of  this  * -neons  1 1  mu  t  ed  (loud 
with  respect,  t.o  aircraft  sampling  time,  it  can  be  seen  (torn  lable 
that  the  calculated  mass  in  .ill  passes  made  prior  *n  puss  1  •,  is  i  nw 
"falling  out"  of  the  cloud,  further,  the  "late  t  ime"  fulluut  (m 
ground)  accounts  for  only  2.2  percent  o<  the  tot  a 1  muss  ln,!cl 

f.) 

(7.9  x  lfV’o)  which  supports  the  observation  that  nr, I,-  .ipn>  i. .  •  ■■  ,r  <■  1  . 


at  Twenty  Minutes  After  Detonation 


15  percent  of  the  cloud's  mass  (as  sampled)  could  be  attributed  to 
particles  >47  pm.  Also  note  that  the  mass  lofted  above  the  visible 
cloud  top  now  accounts  for  only  1  percent  of  the  total  mass  lofted. 

2-4.3  Comparisons  of  the  MBII-1  Dust  Cloud  at  T+10  and  T+20  Minutes 

Figure  2.14  summarizes  the  effects  of  reconstructing  the  MBII-1 
dust  cloud  by  graphically  revealing  how  several  mass  loading  parameters 
are  changing  with  time  and  altitude.  Figure  2.14(a)  depicts  the  dust 
cloud  mass  contribution  from  particles  >47  ,mi  and  shows  where  this  mass 
contribution  is  distributed  at  T+10  and  T+20  minutes  within  the  cloud. 

As  would  be  expected,  much  of  the  mass  contribution  from  these  particles 
in  the  upper  layers  of  the  cloud  at  T+10  minutes  is  being  depleted  with 
time  and  redistributed  throughout  the  lower  regions  of  the  cloud  at 
T+20  minutes. 

Figure  2.14(b)  (plotted  on  same  altitude  scale  as  Figure  2.14(a)) 
depicts  how  the  average  mass  concentration  (all  particle  sizes)  is 
changing  with  time  and  altitude.  From  Tables  2.2  and  2.3  we  can  see 
that  the  mass  aloft  at  T+10  and  T+20  is  nearly  the  same.  However,  due 
to  cloud  diffusion  (cloud  radius  has  increased  by  a  factor  of  2.6), 
Figure  2.14(b)  reveals  an  approximate  factor  of  7  decrease  in  average 
mass  concentration  with  time.  The  similarity  of  the  averaoe  mass  concen¬ 
tration  profiles  at  T+10  and  T+20  minutes,  as  opposed  to  a  "non-linear" 
difference  between  the  curves  in  Figure  2.14(b),  points  out  the  fact 
that  the  MBII-1  dust  cloud  mass  loadino  is  being  governed  by  the  smaller 
particles  aloft,  i.e.,  47  ,.m  pat  *  kies.  Average  mass  concentrations 

from  Table  2.1  (as  sampled)  Pave  been  included  in  Figure  2.14(b)  to 
show  how  'he  cloud  >  ecmist. '  i<  *  int  pi  o<  esses  have  drastically  modified 
the  VP i  *  i  i  a  1  :  I , : j  +  i  1 1 1  i  >  *  * 1  e  "  ea  .u 1  ed  spin.  ’  ra  1  da  *  a  . 


SECTION  3 

THE  MISERS  BLUFF  II  MULTIPLE  BURST  DUST  CLOUD 


3-1  MB  1 1-2  IN-SITU  SAMPLING  DATA  SUMMARY 

Twenty  two  in-si tu  sampling  passes  were  accomplished  by  PMS 
in  the  MBII  multiple  burst  dust  cloud.  Of  these,  sampling  passes  P-11  were 
flown  through  the  cloud  as  the  aircraft  ascended  to  the  cloud  top  and 
are  included  in  the  multiple  burst  data  cloud  analysis. 

Table  3.1  summarizes  and  profiles  the  MBII-2  sampling  mission 
(note  that  the  format  is  identical  to  Table  2.1).  Pass- averaged  mass 
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concentrations  (g/nV )  have  been  calculated  from  the  measured  particle 
number  concentrations  for  each  pass  and  are  tabulated  as  a  function 
of  particle  size,  sampling  altitude,  and  time  of  acquisition.  The  size 
range  of  particles  contributing  to  the  total  mass  concentration  is 
consistent  with  the  size  ranges  used  in  the  MB  1 1 - 1  dust  cloud  analysis 
so  that  mass  loading  comparisons  can  be  readily  made. 

Passes  2,  3  and  4  were  made  through  the  cloud  stem  while  the 
remaining  passes  were  flown  through  the  cloud  cap.  An  interesting 
phenomena  to  note  in  Table  3.1  is  that,  unlike  the  MB  1 1 - 1  cloud  cap, 
mass  concentrations  are  not.  systematically  decreasing  with  time  and 
altitude  but  seem  to  flucuate  randomly  throughout  the  sampling  mission. 
Further,  the  calculated  values  of  mass  in  each  cloud  layer  are  rela¬ 
tively  uniform  throughout  the  cloud  cap  as  opposed  to  those  calculated 
for  the  single  burst  cloud.  The  most  probable  reason  for  these  two 
observations  is  that,  while  under  the  influence  of  oenerally  1  inti t. 
winds  and  little  shear,  wind  induced  diffusion  did  not  play  as  great 
a  role  in  "pulling"  the  MBII-2  dust  cloud  apart  as  was  evidenced  in 
the  single  burst  dust  cloud  (see  section  2-4.3). 

Another  intrinsic  feature  of  the  multiple  burst  dust  cloud 
is  the  larger  percentage  ot  mass  being  contributed  to  t  tie  total  mass 
lofted  by  particles  >4/  .mi  -  approx  ima  I  el  v  P‘>-  IS  percent  in  most  Masses 
as  opposed  to  IS  percent  in  the  single  hurst  cloud.  Ihoro  are  prohaMv 


Table  3.1  MISERS  BLUFF  1 1 - 2  In-Situ  Sampling  Data  Summary 


two  contributing  factors  influencing  this  observation:  (1)  MB  1 1 -2  flow 
fields  were  sufficient  to  keep  these  larger  particles  aloft  for  sustained 
periods  of  time  and;  (2)  the  MBI 1-2  ANFO  stack  site  geologies  varied 
somewhat  from  the  single  burst  GZ  site  geology. 

The  MBI 1-2  dust  cloud  mass  loading  analysis,  based  on  the 
same  assumptions  made  for  the  single  burst  cloud  analysis  (Section  2-1), 

revealed  that  the  total  mass  lofted  in  the  multiple  burst  dust  cloud 
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was  5.1  x  10  grams.  Again,  this  value  compares  favorably  with  5.2  x  10' 

grams  determined  by  PMS  in  their  independent  analysis. 

3-2  MBI 1-2  DETAILED  SAMPLING  PASS  DATA 

Figures  3.1-3.11  in  this  section  are  identical  in  format  to 
those  presented  in  Section  2-2.  For  each  sampling  pass  used  in  the 
analysis,  pass- averaged  particle  size  and  cumulative  mass  distributions 
have  been  plotted  (Figures  3 . X ( a ) )  along  with  calculated  mass  concen¬ 
trations  within  the  cloud  (Figures  3 . X ( b ));  except  sampling  Pass  1  which 
was  flown  below  the  cloud  cap  and  outside  the  visible  stem.  In  the  particle 
size  distributions,  circles  represent  particles  •  47  pin,  squares 
represent  particles  47-470  pm,  and  triangles  depict  the  size  distri¬ 
bution  of  particles  ■  470  pin.  For  the  mass  concentration  plots,  the 
contribution  from  particles  47  pin  are  plotted  with  circles  (solid 
line)  and  that  from  particles  47  are  plotted  with  squares  (dashed 
line). 

3-3  GENERAL  COMMENTS  CONCERNING  THE  MBI 1-2  DUST  CLOUD 

3-3.1  Particle  Size  Distributions 

As  in  the  MB  1 1  - 1  dust  cloud,  particle  sizes  47  ,.m  are 
approximately  distributed  accord  inn  to  a  power  law  (see  Section  2-3.1), 
however,  there  appears  to  be  more  variation  of  the  power  "q"  in  the 
multiple  burst  dust  cloud.  For  the  part  icle  size  ranee  47-470  .  m, 

3.0  <1  5.0 

while  for  particles  470  . m  (except  passes  3  and  9), 
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q  5.0 


Again,  it  should  be  mentioned  that  the  size  distribution  depicted  in 
Figures  3.X(a)  are  based  on  pass- averaged  data  and  significant  varia¬ 
tion  in  "q"  will  exist  from  point  to  point  within  a  given  pass  and 
in  the  tail  end  of  the  distributions  due  to  poor  counting  statistics. 

3-3.2  Cumulative  Mass  Distributions 

Passes  2  and  3  were  flown  through  the  visible  stem  of  the  cloud  and 
the  cumulative  mass  distributions  show  that  the  dust  mass  is  fairly 
evenly  distributed  between  particles  less  than  and  greater  than  47  uni. 

At  later  sampling  times,  the  major  contribution  to  mass  loading  shifts 
toward  the  smaller  particles  where  approximately  65-75  percent  of  the 
mass  is  found  (except  passes  9  and  11  where  the  contribution  from 
particles  <47  pm  is  80-85  percent). 

3-3.3  Mass  Concentration  Time  Histories 

Unlike  the  MBII-1  cloud  penetrations,  the  highest  mass  con¬ 
centration  encountered  in  the  MB  1 1 -2  cloud  was  in  the  cloud  cap 
(pass  6)  and  not  in  the  cloud  stem.  One  reason  for  this  is  that,  because 
of  operational  problems,  cloud  penetrations  in  the  single  burst,  cloud 
began  approximately  two  minutes  before  entry  into  the  multiple  burst 
cloud  was  possible,  i.e.,  when  stem  mass  concentrations  were  much 
higher.  Another  possibility  is  that  flow  fields  within  the  multiple 
burst  dust  cloud  were  Sufficient  to  keep  those  portions  of  the  cloud 
containing  the  higher  mass  concentrations  aloft,  for  a  sustained  pet  iod 
of  time,  whereas,  the  mass  in  the  MBII-1  dust,  cloud  was  rapidly  bo  inn 
depleted  by  fallout  through  the  stem  region  at  very  early  times.  I'oal 
mass  concentrations  in  the  MB  1 1-2  cloud  were  generally  3-4  times  the 
average  values;  pass  5  had  a  peak  value  slightly  5  times  the  avetage 
mass  concentration. 


Figure  3.1.  MISERS  BLUFF  1 1  - 2 ,  Pass  1,  Mass  Concentration  Time  History 
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Figure  1.2(a).  MIS F RS  BLUFF  11-2,  Pass  2,  Particle  Size  and 
Cumulative  Mass  Distributions 
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.  MISERS  BLUFF  1 1 -2 ,  Pass  4,  Mass  Concentration  Time  History 
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Particle  Diameter  (;;m) 

Figure  3.5(a).  MISERS  BLUFF  II-Z,  Pass  5,  Particle  Size  and 
Cumulative  Mass  Distributions 


Number  Concentration  (Particles/cm 


Pass  5,  Mass  Concentration  Time  History 
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re  3.6(b).  MISERS  BLUFF  II-2,  Pass  6,  Mass  Concentration  Time  History 
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RECONSTRUCTING  THE  MISERS  BLUFF  1 1-2  DUST  CLOUD 


In  a  manner  identical  to  that  used  in  the  MB 1 1 - 1  dust  cloud 
reconstructions  (see  Section  2-4  for  methodology  and  assumptions), 
the  multiple  burst  dust  cloud  was  reconstructed  at  ten  and  twenty 
minutes  after  detonation  to  obtain  "snapshots"  of  the  entire  cloud 
at  several  times  for  comparison  purposes. 

3-4.1  The  MBII-2  Dust  Cloud  at  Ten  Minutes  After  Detonation 

Table  3.2  is  the  summary  of  calculated  mass  values  within  the 
T+10  minute  dust  cloud  as  a  function  of  altitude.  Average  mass 
concentrations  have  been  calculated  by  summing  mass  contributions  from 
all  sized  particles  at  a  given  altitude  and  then  dividing  by  the 
estimated  cloud  volume  at  that  altitude.  Values  of  cumulative  mass 
have  been  tabulated  to  insure  that  all  cloud  mass  has  been  conserved 
in  the  reconstruction  process. 

As  in  the  single  burst  dust  cloud  at  T  + 1 0  minutes,  some  of  the 
clouds  mass  (  7  percent)  appears  above  the  top  of  the  visible  cloud  (see 
Section  2-4.1  for  a  discussion  on  the  possible  explanation  for  this  con¬ 
tradiction  between  the  inferred  T'lf)  minute  cloud  and  the  visible  (lend: 
Since  this  "phenomenon"  is  common  to  both  dust  clouds  at  T-*l'i  minute-,,  1 
is  not  felt  that  the  mass  loading  ratios  would  be  s  i  on  i  t  i  can '  1  v  dHn'ol 
by  it. 

In  order  to  smooth  out  some  of  the  fluctuation  in  aveiaoe  mans 
concentrations  in  the  cloud  cap,  passes  /  and  K  and  passes  u  and  1"  nave 
been  averaged.  The  calculated  averages  appear  on  the  <  ight-hand  side  ni 
Table  3.2.  Similarly,  pass  data  from  passes  ?,  3,  and  4  nave  beet 
averaged  in  calculating  an  average  mass  concent  i  at  ion  in  the  >.  load  s'e<". 

3-4.2.  The  MBII-2  Dust  Cloud  at  Twenty  Minutes  After  lie t nna ' i on 

Table  3.3  sunimar  i  /(.*•»  t  no  gross  tea?  sees  of  *  he  mill  *  1  ?•  le  Lea' 
dust  cloud  at  T'20  minutes.  2/  this  ?  ime,  or  I v  2.  pe> <  ei  ?  of 
mass  lotted  appears  above  t  tie  visible  (loud  'op.  I  w  "la'e 

r  i  i  Me"  fa  limit,  i  .  e  .  ,  *  he  n  ass  t  ha  *  is  ot  t  n*.  mound  •  i  Table  ■ .  , 


Table  3.2  MISERS  BLUFF  1 1 -2  Dust  Cloud  Reconstructed  at  Ten  Minutes  After  Detonation 


Table  5.3  MISERS  BLUFF  1 1 -2  Dust  Cloud  Reconstructed  at  Twenty  Minutes  After  Detonation 


accounts  for  r6  percent  of  the  mass  lofted.  Recalling  that  at  T+20 
minutes,  the  "late  time"  fallout  from  the  MBII-1  cloud  amounted  to 
only  2.2  percent  of  the  lofted  mass,  the  6  percent  value  for  MB  1 1  - 2 
supports  the  observation  that  a  much  greater  percentage  of  the  cloud's 
mass  was  in  the  form  of  particles  that  could  "fall  out"  of  the  cloud, 
i .e. ,  particles  y  47  pm. 

3.4.3  Comparisons  of  the  MB  1 1-2  Dust  Cloud  at  T+10  and  T+20  Minutes 

Figure  3.12(a)  shows  how  particles  47  ,.m  are  being  redis¬ 
tributed  in  the  MB  1 1 - 2  dust  cloud  as  a  function  of  altitude  and  time- 
after-burst.  As  in  the  MBII-1  cloud,  the  upper  cloud  layers  are  being 
depleted  with  time  with  a  corresponding  mass  enhancement  in  the  lower 
regions  of  the  cloud. 

Average  mass  concentrations  within  the  MB  1 1-2  dust  cloud  at 
T+10  and  T+20  minutes  have  been  plotted  in  Figure  3.12(b)  as  a  function 
of  altitude.  Total  mass  aloft  at  T+20  minutes  is  nearly  the  same  as 
that  aloft  at  T+10  minutes  (see  Tables  3.2  and  3.3),  while  Figure  3.12(b) 
shows  a  factor  of  2-3  decrease  in  average  mass  concentration  with  time. 
(Ten  minutes  of  cloud  diffusion  resulted  in  a  factor  of  1.6  increase 
in  the  cloud  radius  at  T+20  minutes;  compare  this  with  the  2.6  increase 
in  cloud  radius  for  the-  highly  sheared  MBII-1  dust  cloud.)  'Ji  thin  the 
same  time  frame,  there  was  a  factor  of  7  decrease  in  the  average  mass 
concentration  in  the  MBII-1  dust  cloud  (Section  2-4.3).  Tin's  further 
supports  the  observation  that  the  cingle  burst  dust,  cloud  was  strongly 
affected  by  significant  variations  in  wind  speed  and  direction  aloft, 
(shear)  that  we-'e  not  evident  in  the  wind  profile  obtained  during  the 
multiple  burs,  event. 
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SECTION  4 

SINGLE  AND  MULTIPLE  BURST  DUST  CLOUD  COMPARISONS 
4-1  INTRODUCTION 

In  current  multiple-burst  dust  cloud  modeling,  the  total  mass 
aloft  at  a  given  time  after  stabilization  is  assumed  to  be  a  simple 
addition  of  the  masses  calculated  for  the  corresponding  single-burst 
dust  clouds.  Therefore,  if  a  dust  cloud  resulting  from  a  single 
detonation  contains  X  mass,  the  stabilized  dust  clnuH  from  M  identical 
detonations  will  contain  NX  mass  at  the  same  time  after  detonation. 

The  stabilizeu  (T+10  minute)  clouds  for  MISERS  BLUFF  II  were 
loaded  almost  exclusively  by  very  small  particles.  These  high  drag 
particles  will  tend  to  be  contained  in  the  visible  cloud  whose  volume 
is  dominated  by  the  cloud  cap.  To  take  advantage  of  this  result,  only 
the  MBII  cloud  caps  will  be  compared  at  similar  times  after  detonation. 

4-2  MBII  MASS  LOADING  RATIOS 

Because  most  of  the  assumptions  that  were  made  in  calculating 
total  lofted  mass  (e.g.,  dust  specific  gravity,  spherical  particles, 
and  circular  cloud  geometry)  were  the  same  for  both  the  single  and 
multiple  burst  dust  clouds,  multiplicative  uncertainties  attributable 
to  these  assumptions  should  cancel  out,  or  tend  to  become  less 
important,  when  mass  loading  ratios  are  calculated.  These  and  other 
uncertainties  are  discussed  in  more  detail  in  Section  5. 

From  Tables  2.2  and  3.2  (MBII  dust  clouds  at  T+10  minutes), 
calculated  values  of  cumulative  mass  above  the  cloud  stem  show  that 
0.64  x  10^  kg  and  4.5  x  10^  kg  of  dust  remains  aloft  in  the  sinole 
and  multiple  burst  cloud  caps,  respectively,  at  cloud  stabilization. 

The  mass  loading  ratio  between  the  two  dust  clouds  at  this  time  is 
approximately  7.1.  By  twenty  minutes  after  detonation,  the  single  and 
multiple  burst  dust  clouds  contain  0.63  x  10^  kg  and  4.15  x  h/’  kg  of 
dust,  respectively,  yielding  a  mass  loading  ratio  of  6.6.  The  downward 
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trend  of  the  mass  loading  ratio  with  time  is  due  to  the  fact  that  the 
mass  of  the  multiple  burst  cloud  is  comprised  of  a  greater  percentage 
of  particles  -47  pm  which  are  "falling  out"  of  the  cloud  cap.  If  a 
similar  analysis  were  to  be  extended  out  to  times  greater  than 
twenty  minutes,  it  would  be  apparent  that  the  mass  loading  ratio 
between  the  cloud  caps  asymtotical ly  approaches  a  value  of  6.0,  i.e., 
the  ratio  of  masses  in  the  cloud  caps  being  determined  only  by  particles 
•47  pm.  Within  the  range  of  measurement  uncertainties,  these 
calculated  values  of  mass  loading  ratio  ajapear  to  approximate  what 
would  be  predicted  by  linear  superposition,  i.e.,  six. 

Recalling  that  the  T+10  minute  dust  cloud’s  mass  loading 
ratio  was  calculated  to  be  >-7.1,  measurement  uncertainties  aside, 
there  is  some  evidence  that  the  18  percent  "enhancement"  (7.1/6)  in 
mass  loading  could  be  real.  First,  the  calculated  enhancement  in 
the  mass  loading  ratio  is  governed  by  a  factor  of  12  increase  in 
mass  aloft  of  particles  47  pm  in  the  M B 1 1-2  dust  cloud;  a  factor 
of  six  would  have  implied  linear  superposition.  While  differences 
in  site  geology  (charges  3  and  4  were  detonated  over  a  silty/sandy 
surface  layer  approximately  0.5m  thick;  see  Ref.  6)  could  account 
for  this  phenomenon,  photographic  records  of  the  multiple  burst  event 
(Ref.  5)  strongly  suggest  that  an  enhanced  flow  field  did  exist  that 
could  be  responsible  for  lofting  the  large  particles  in  the  multiple 
burst  cloud.  Hence  there  was  not  only  a  source  of  large  particles 
present  for  the  MB  1 1  - 2  dust  cloud,  but  there  also  appears  to  be 
a  mechanism  present  to  loft  these  larger  particles  not  evident  in 
the  MB  1 1 - 1  cloud  photography  -  enhanced  flow  fields. 

Reiterating,  the  18  percent  increase  in  mass  loading  ratio 
could  be  attributable  to  measurement  uncertainties,  howevei .  there 
is  st  rong  evidence  that  enhanced  flow  fields  generated  by  i ntei act  i  no 
shockwaves  may  also  be  a  factor. 
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MB 1 1  CLOUD  VOLUME  RATIO 


In  addition  to  predicting  a  six-fold  increase  in  mass  aloft 
at  cloud  stabilization,  linear  superposition  implies  that  dust  mass 
concentrations  (g/m  )  are  additive  in  regions  where  the  dust  clouds 
overlap.  For  instance,  if  six  detonations  occurred  simultaneously 
at  the  same  location,  linear  superposition  would  predict  a  multiple 
burst  dust  cloud  that  had  the  same  dimensions  as  a  single  burst  dust 
cloud  but  would  contain  six  times  the  mass.  This  would  imply  that  the 
average  mass  concentration  within  the  stabilized  multiple  burst  dust 
cloud  would  be  six  times  that  found  in  the  single  burst  cloud  at  the 
same  time-after-burst. 

Figure  4.1  shows  how  some  measured  dimensions  of  the  MBII 
multiple  burst  dust  cloud  (dashed  lines)  compare  with  those  cloud 
dimensions  that  would  be  predicted  by  linearly  superimposing  six  of 
the  MBII  single  burst  clouds  (solid  lines)  at  ten  minutes  after 
detonation.  Figure  4.1(a)  depicts  the  vertical  cross-section  of  the 
measured  and  predicted  multiple  burst  dust  clouds  while  Figure  4.1(b) 
is  an  overhead  (planar) view  of  the  two  dust  clouds.  The  hexagonal 
array  of  dots  in  Figure  4.1(b)  is  the  multiple  burst  charge 
configuration  and  has  been  drawn  to  scale  (200  meter  centerline 
di ameter ) . 

One  immediate  observation  is  that  the  measured  volume  of  the 
multiple  burst  cloud  is  significantly  greater  than  would  be  predicted 
by  linear  superposition.  Calculations  based  on  cloud  diameter  and 
thicknesses  at  T+10  minutes  (Tables  2.2  and  3.2)  reveal  that  the 
measured  multiple  burst  cloud  volume  is  2.7  times  that  predicted 
by  superimposing  six  of  the  single  burst  clouds.  Assuming  that,  the 
predicted  multiple  burst  cloud  contains  3..''.  x  10fl  kg  of  dust  (six 
times  the  mass  in  the  single  hurst  cloud),  linear  superposition 
would  predict  an  average  mass  concentration  in  the  multiple  burs* 
dust  cloud  of  0.67  o/m’.  The  average  mass  concent i at  inn  in  the 
multiple  burst  dust  cloud  was  measured  as  on  1  y  f  i .  •:  o'n  'sen 
Table  3.2). 


Another  interesting  parameter  to  calculate  is  the  total  mass 

that  would  probably  be  encountered  by  a  missle  ascending  (or  descending) 

vertically  through  the  stabilized  dust  clouds,  i.e.,  / pdh .  Assuminci 
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that  a  cross  sectional  area  of  1  m  is  being  swept  by  the  "missile", 
calculations  show  that  linear  superposition  would  predict  a  value  of 
integrated  mass  in  the  six  MB  1 1  single  burst  clouds  to  be  approximately 
twice  that  calculated  for  the  MB  1 1  multiple  burst  dust  cloud  (  1340  a 
versus  720  g).  It  appears  then  that  linear  superposition  of  dust 
clouds,  formed  from  closely  spaced  multiple  detonations,  tends  to 
overestimate  the  dust  environment  (g/m  ).  Further,  this  overestimate 
in  mass  concentration  is  due  mainly  to  a  corresponding  underestimate 
in  the  multi-burst  cloud  volume  predicted  by  linear  superposition. 
However,  if  the  number  of  burst  in  a  given  attack  scenario  is  very 
large  and  closely  spaced, and  they  encompass  a  relatively  large  area 
(unlike  MB  1 1  -  2 ) ,  the  central  clouds  will  be  radially  confined  and 
mass  concentrations  within  this  area  would  not  drop  off  as  sharply  as 
indicated  by  MB  1 1-2  (approximately  a  factor  of  2).  This  would  represent 
a  "worst  case"  situation  and  linear-  superposition  would  probably 
yield  very  good  results. 

In  addition  to  the  apparent  weaknesses  of  linear  superposition 
in  describing  the  multiple  burst  environments  that  have  already  been 
mentioned  in  this  section.  Figure  1.1(a)  reveals  another  potential 
problem  area  --  that  of  inadequately  describing  the  altitude  and 
thickness  of  the  multiple  burst  dust  cloud.  This  aspect  has  also 
been  raised  try  unpublished  hydrocode  (HULL)  calculations.  The 
primary  cause  is  an  energy  superposition,  i.e.,  six  times  the  ener gy 
went,  into  only  7  times  the  volume;  the  excess  energy  density  was 
probably  used  primarily  in  raising  the  multiple  burst  (.loud  to  a 
higher  alti'ude.  It  j  •,  not  felt  that  Hie  slight  ditteteiK.es  in 
atmospheric  temper  a t u re  lapse  rates  existing  at  the  turn*  of  each 
detonation  contributed  to  this  difference  in  (.loud  altitudes  --  in 
frtc*.  the  air  was  more  stable  din  ino  the  multiple  burst  test. 


SECTION  5 

UNCERTAINTIES  IN  THE  MBII  DUST  CLOUD  MASS  CALCULATIONS 


5-1  INTRODUCTION 

The  need  to  discuss  uncertainties  became  particularly  obvious 
when  it  is  realized  that  the  calculated  values  of  total  lofted  mass 
during  each  of  the  two  events  in  the  MISERS  BLUFF  II  program,  exceeded 
a  nuclear  cloud  "hydrodynamic  limit"  prediction.  This  "limit"  is 
based  on  hydrocode  flow  field  calculations  (SHELL)  which,  for  para¬ 
metric  nuclear  dust  cloud  simulations,  limits  the  total  lofted  mass  to 
approximately  one  ton  of  dust  per  ton  of  nuclear  yield.  Due  to  the 
way  net  energy  output  is  partitioned  in  nuclear  devices  (i.e.,  much 
of  the  energy  being  expended  in  the  form  of  radiant  energy),  it  is 
common  practice  to  equate  1  ton  of  TNT  yield  to  2  tons  of  nuclear  yield 
for  blast  simulations.  The  hydrodynamic  limit  for  high  explosive  (HE) 
detonations  is,  therefore  (for  an  exactly  analogous  flow  field),  two 
tons  of  dust  per  ton  of  TNT  yield. 

In  the  MBII  single  burst  detonation  (10(1  tons  TNT  equivalent), 
the  total  lofted  mass  was  calculated  to  be  7.9x10  kg  (Section  2-1)  or 
B.B  tons  of  dust  per  ton  of  yield.  Similarly,  total  lofted  mass  in 
the  multiple  burst  event  (600  tons  TNT  equivalent)  was  calculated  to 
be  5.1x1 0^*  kg  (Section  3-1)  or  9.4  tons  of  dust  per  ton  of  yield  — 
clearly  much  more  than  "allowed"  bv  'he  HI  hydrndvrnmii  limi'.  Ihe  im¬ 
pose  of  this  section  is  to  discuss  some  of  the  assumptions  that  wen  e 
made  in  calculating  these  values  of  MBII  lofted  mass  and  assess  t he* i * 
effect  on  tin.'  calculations.  Further,  the  validity  of  using  the  pin  Mm* 
hydrodynamic  limit  to-  limiting  dust  lofted  in  HI  detonations  is 
discussed  in  Section  5-3. 

5-2  SOllP.Lf  5  IMP  I  INF!  P  TAINT  If  S 

5-2.1  In-Si’u  Sampling  IHofiles 

Ihe  ii  i  i  bn  i  r  ,e  dus '  (loud  sat'i'i  1  i  nq  missions  we  •  <■  <  1  '■  g'  '  planed  ' 

'  t  ” ' ; ;  If'  !)(.•(■  - .  ]  f  i  i  j .  i  • ,  t  >  i  r  !  ii  ■  1  <  •  *.<•  p.r.e  'is  fluid  *  i .;  •  •  s  (-.'<••  -  i  i  i  i  .i  s  ;  i  ; 

sampling  pass  al’i-udes.  When  pfcssibln,  sampling  passes  we*  <•  cade  u 


generally  up-wind/down-wind  directions  so  that  the  aircraft  traversed 
regions  of  subvisible  fallout  from  higher  altitudes  for  sheared  clouds. 
The  main  purpose  for  having  identical  flight  plans  for  both  sampling 
missions  was  to  insure  a  common  data  base  for  subsequent  dust  cloud 
reconstruction  efforts  and  mass  loading  comparisons.  Ttie  major  draw¬ 
back  in  the  flight  plans  was  that  by  sampling  upwind/downwind,  the 
largest  "dimension"  of  the  dust  clouds  was  sampled  at  each  pass  altitude. 

One  of  the  assumptions  made  in  the  mass  calculations  was  that 
the  dust  cloud  was  circular  at  each  sampling  altitude  and  t. fie  sampling 
pass  length  was  the  diameter  of  tins  circle.  Clearly  this  assumption 
could  lead  to  overestimating  the  mass  at  altitude,  especially  it  the 
clouds  had  become  elliptical  due  to  wind  induced  shear  and  diffusion. 
Since  all  of  the  sampling  passes  were  not  made  upwind/downwind  (Ref.  2), 
a  significant  effort  was  expended  in  analyzing  the  sampling  path  length 
data  and  recomputing  the  cloud  volume  at  each  altitude  assuming  the 
cloud  cross  section  was  elliptical  rather  than  circular.  Tins  effort 
resulted  in  decreasing  the  total  mass  lofted  in  the  single  and  multiple 
burst  clouds  by  21  percent  and  17  percent,  respectively.  To  he  con¬ 
servative,  let  25  percent  be  an  upper  limit,  for  the  mass  overestimate 
due  t.o  assumptions  concerning  cloud  geometry. 

Another  aspect,  of  ttie  in-situ  sampling  profiles,  which  could 
lead  to  significant  errors  in  ostimatino  the  total  mass  lofted,  is  t he 
time  after  detonation  that  initial  penetrations  of  the  clouds  wore  made, 
lor  reasons  net  •  a  i  n  i  i  o  >o  c  >  ew  /  f<>'  he*  (light  npe>  a  r  ions  on  uomc', 

in-si ‘u  samel  i  no  began  much  earlier  in  'tie  single  hurst  <  loud,  when  "ass 
v.r.r  ■  I'fit.i  at  ii'ir.  wi""  V*  mv  hi  :li ,  than  w.r.  !><■•.•.  i!>  !<•  in  the  -  ulf  i  j  1 1  *  >  !<',n  •  t 
■  !  i  cl  I .  eon  o  .  /  tin  ap"r'i  nil  n  t  .it  "h  in  t  iorr.  |  .mi  |  -  I  ,  • 

I  ' '  t  f  > "  f  in  Cue  i  j  I  t  i  n  |  e  burnt  !  >  :i  i  m  I  ‘  1  <  n  I  •  I  •  ■  <  e  .  t  e  i  1 1  !  i '  n  ‘  i .  ea  •"  1  . 
wS-en  Ml  Ion*. 
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Dust  Particle  Shape  and  Specific  Gravity 

Throughout  the  entire  comparati ve  analysis,  the  dust  particles 

have  been  assumed  to  be  spherical  in  shape  and  have  a  specific  density 
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of  2.3  g/cm  .  Significant  overestimates  in  calculated  values  of  lofted 
mass  could  be  incurred  if  the  dust  particles  were,  in  fact,  a  totally 
different  shape.  Since  most  of  the  mass  lofted  is  attributable  to 
particles  -47  pin  in  diameter  (see  Figure  4.1),  this  discussion  will 
focus  on  particles  in  that  range  of  tfie  size  spectrum. 

At  DNA's  request  for-  a  secondary  mission,  the  Air  Torce  Technical 
Applications  Center  (AFTAC)  flew  three  sampling  passes  through  the  top 
of  the  MB  1 1  multiple  burst  dust  cloud  between  fifteen  and  twenty  minutes 
after  detonation.  During  each  sampling  pass,  AFTAC  obtained  filter 
paper  samples  containing  dust  from  the  cloud  and  one-half  of  each  filter 
was  subsequently  forwarded  to  CAI  for  analysis.  The  particular  analysis 
performed  on  the  samples  was  called  a  Suspended  Particulate  Evaluation 
and  Classification  (SPEC)  (it  was  accomplished  by  Materials  Consultants 
&  Laboratories,  Inc.  of  Monroeville,  PA).  The  output  from  the  SPEC 
analysis  is  a  detailed  description  of  particle  size,  shape,  and 
chemistry  (elements)  for  the  approximately  1000  particles  in  em  h 
'•ample  analyzed.  With  respect  to  the  particle  shape  of  dust  particles 
1  .  the  analysis  revealed  a  most  probable  wi dth-to- length  ratio 
of  rl-r>M,  i.e.  ,  a  lenqt.h-tn-width  ratio  (r)  of  1.7.  (Moan  0.5b  and 
Standard  deviation  0.17.)  Assuming  the  "real"  shape  of  the  nartii  les 
wen'  i  irculur  discs  (see  Append i«  0)  with  an  aspect  ratio  (r)  of 
l./,  '"ass  Tabulations  that  were  'ade  assuming  Particle  sphermity 
would  overest.  j- at>>  the  mass  lofted  h/  Is  percent  (4  r/f. ) . 

1  he  el  i  At  filter  samples  also  provided  ',A !  with  a  totally 
independent-  •  ass  •  oasureneri  t.  free  the  Ml;  II -  2  duct  t  In  el.  by 

1  'P'W  ill  O  file  1  y  ‘ 1  n  I  i  e  Weight  of  the  (life,  ;  ape  >  -  |  a  j  , ,  to  i  b  ,  S'  I 
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the  pass  averaged  mass  concentration.  AFTAC  furnished  SAI  with  their 
best  estimate  of  the  filter  paper  tare  weights  and  individual  pass 
Sample  volumes.  The  "dirty"  filter  papers  were  weighed  by  SAI  and 
mass  concentrations  were  calculated  for  each  sampling  pass.  Table 
5.1  is  a  tabulation  of  this  data. 

Table  5.1  MB  11-2  Dust  Cloud  Hass  Concentrations  Inferred 
From  AFTAC  Filter  Paper  Samples 

Start  Altitude  Mass 

Pass  No.  Time  (MST)  MSL  (m)  Concentration  (g/in  ) 

1  13:14:57  3540  .094  *.02 

2  13:17:16  3410  * 

3  13:20:10  3350  .11  *  .04 

*The  dust  contribution  to  the  filter  paper  weight  was  not  discernable 
using  the  best  estimated  tare  weight.  Filter'  paper  appeared  much 
"cleaner"  than  the  others;  probably  penetrated  patchy  portions  of 
the  cloud  cap. 

With  the  exception  of  sampling  pass  2,  the  AFTAC  data  i nd ■>  .ates  an 
average  dust  concentration  within  tire  top  of  the  MBII-2  oust  cloud  of 
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approximately  0.1  g/m  at  15-20  minutes  after  detonation.  Comparing 
this  number  with  that  calculated  for  the  same  altitude  layer  (3140- 
3400m)  iri  the  reconstructed  MB  1 1 -2  dust  cloud  (Table  3.3  or  figure  ml 
we  see  excellent  agreement. 

The  SPEC  analysis  of  particle  sizes  ■  4/.  m  rrom  ’he  AFTAr 
filters  also  revealed  that  30  percent  of  tin'  total  sample  mass  was 
contributed  by  particles  ‘  15  ,.m.  This  is  in  very  good  agreement,  with 
t  fie  special  da!  a  obtained  bv  PM5. 

Again  it  is  encouraging  to  find  such  good  agreement,  in  the 
results  obtained  from  two  totally  independent,  data  bases .  based  mi 
the  SPL(.  ass  analysis  and  the  cnloulatd  values  of  mass  <  one  mil t  af  ion 
fror  the  Af  T  Af"  filter  napers,  it  is  felt  that  the  assir  ;d  i  ore,  made 
concet  n  i  n”  ’he  dust  particles  specif  it  gravity  being  w<r  a  'easona 
ore  and  did  no*  lead  In  n  vet  eg  t  i  ma  t  i  ng  ti,e  *o’al  irass  letted. 
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5.3 


VALIDITY  OF  THE  NUCLEAR  HYDRODYNAMIC  LIMIT  FOR  HE  TESTS 

Based  on  the  aforegoing  discussion  concerning  possible  sources 

of  uncertainties  in  the  MB  1 1  dust  cloud  mass  calculations ,  it  seems 

★ 

one  could  justify  an  approximate  40  percent  reduction  in  the  calculated 
values  of  lofted  mass  for  each  event.  This  reduction  would  result  in 
a  lofting  ability  of  5.3  and  5.6  tons  of  dust  per  ton  of  TNT  yield  for 
the  single  and  multiple  burst  event,  respectively  -  still  much 
greater  than  the  2  tons  "permitted"  by  the  HE  hydrodynamic  limit. 

Exceeding  the  hydrodynamic  limit  is  not  unique  to  the 
MB  1 1  tests.  For  HE/ANFO  event  DICE  THROW  (500  tons  TNT  equivalent),  it 
has  been  estimated  that  6.0  tons  of  dust  were  lofted  (above  the 
initial  fireball  top)  per  ton  of  TNT  yield  (Ref  4),  hence,  the  nuclear 
hydrodynamic  limit  may  not  be  applicable  to  HE/ANFO  dust  clouds.  One 
possible  reason  for  this  is  that  the  dust  loading  mechanisms  for 
nuclear  and  HE  fireballs  are  very  different.  Nuclear  fireballs  are 
loaded  with  dust  preferentially  near  the  center  line  of  the  fireball 
flow  field  -  in  fact  the  parametric  SHELL  calculations  commenced  witti 
all  mass  being  initially  placed  in  the  inner  one-eioht h+of  the  hemis¬ 
pherical  fireball  volume  (surface  burst).  At  this  location  initially, 
much  of  the  mass  would  be  lost  immediately  due  to  the  base  surer 
tendency  along  the  centerline  and  would  not  be  "available"  tor 
subsequent  lofting.  An  HE  firehill  (ANI'O-i  upped  <ylindei  ihun'ie), 
however,  is  loaded  by  eieit.u  preferentially  outside  the  .  entei  of 
the  fjrehtil  1.  being  loaded  in  this  region,  the  loss  o<  mass  ulor'i 
trie  stagnant  renter  axis,  as  ptediefed  hy  'Ur  puruneft  i<  ‘•HI  I  ! 
i.  i  lull  at  ions,  does  not  m.t  lit  and  more  ass  is  "ova  i  1  at  |e '  to  be 
lofted.  t  P(M  all  fh.lt.  ft  o  III  I  a  I  e  !  :U  >  t  ,  'finny  all  t  flow  is  i  n  i  t  i  a  f  "C 
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predicted  by  SHELL  calculations  for  a  nuclear  dust  cloud.  While  the 
buoyant  energy  is  comparable,  the  hydrodynamic  flow  fields  of  both 
the  ejecta  and  the  fireball  are  very  different  for  nuclear  and  HE/ANrO 
bursts . 


SECTION  6 

CONCLUSIONS  AND  RECOMMENDATIONS 


An  indepth  analysis  of  the  data  obtained  during  Phase  1  of 
the  MISERS  BLUFF  II  Cloud  Sampling  Program  indicates  that  the  MB  1 1 
multiple  burst  dust  cloud  could  not  be  completely  modeled  by  linearly 
superimposing  six  of  the  MBII  single  burst  clouds.  As  might  be 
expected,  the  major  weakness  in  the  linear  superposition  modeling 
approach  is  its  inability  to  treat  the  effects  of  the  hydrodynamic  and 
thermodynamic  interactions  occurring  between  the  shock  fronts/fireballs 
of  a  limited  number  of  closely  spaced,  simultaneous  detonations. 

As  a  result  of  this  weakness,  several  conclusions  can  be  made 
based  on  the  analysis  performed  in  this  report  for  the  MBII  dust 
clouds : 

•  Linearly  superimposing  six  of  the  MBII  single  burst  dust, 
clouds  at  T+20  minutes  resulted  in  oy_e_re_s_tim_a_t i jkj  the 
vertically  integrated  dust  mass  aloft  (,'rdh)  by  a  fact.oi 
of  2  over  the  measured,  multi  burst  dust  environment. 

This  result  is  due  mainly  to  the  inability  of  linear  super¬ 
position  to  accurately  predict  the  volume  of  the  dust  cloud 
formed  from  a  few  closely  spaceu,  interacting  detonations. 

a  While  overestimating  the  magnitude  of  the  multiburst  dust 
environment  (mass  concentrat ions ) ,  linear  superposition 
appears  to  underestimate  the  areal  extent  of  the  actual 
environment  that  is  produced  in  a  closely  spaced,  multi¬ 
burst  situation.  Again,  this  conclusion  is  a  result,  of 
linear  superposition  not  treating  thermal  energy  diss ina ' ion 
and,  hence ,  incorrectly  predict  inn  mul  t  iln*  <  Lnel  vrl  nr  os. 

A  Within  the  range  of  measurement  uncertain;  it's  ,  linen*  sunn  - 
position  may  nredic t  reasonable  values  of  mass  nlnt‘,  i.e., 
from  N  identical  simultaneous  detonations  (In'  cm  1 
loft  <  mass .  However,  toe  tart  1 1  .at  a  "se  m*r  ■ 
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centra]  plume  of  dust  was  evident  in  the  MBII-2  cloud 
photography,  first  visible  approximately  20  seconds  after 
detonation,  provides  strong  evidence  that  the  measured 
^7X  increase  in  mass  aloft  in  the  multiple  burst  cloud 
may  be  very  real . 

Based  on  the  above  conclusions  concerning  the  validity  of 
utilizing  linear  superposition  in  defining  dust  clouds  from  a  few 
closely  spaced,  multiple  burst  detonations,  the  following  recommen¬ 
dations  are  made: 

•  Cloud  volume  calculations  were,  for  the  most  part,  derived 
from  aircraft  sampling  pass  lengths  at  or  very  near  the 
times  of  dust  cloud  reconstructions  (T+10  and  T+20  minutes). 
More  accurate  values  could  be  obtained  by  a  detailed 
analysis  of  all  the  MBII  dust  cloud  photography  as  taken 
from  several  azimuths.  Because  cloud  volume  calculations 
were  instrumental  in  arriving  at  several  of  the  conclusions, 
it  is  recommended  that  aJJ_  photographic  documentation 

(TIC  and  SRI)  of  the  MBII  dust  clouds  at  T+10  minutes  be 
analyzed  to  derive  more  accurate  values  of  cloud  volume. 

•  Some  future  analytical  effort  should  be  expended  to  under¬ 
stand  how  interacting  fireballs  dissipate  their  buoyant 
energy  as  they  rise  in  the  atmosphere.  This  information 
is  needed  to  accurately  predict  dust  cloud  volumes  and, 
hence,  mass  concentrations  within  any  dust  cloud  generated 
by  closely  spaced,  simultaneous  detonations  where  linear 
superposition  of  cloud  volume  is  not  valid.  The  buoyant 
energy  partitioning  (for  a  suddenly  released  volume  of 
buoyant  fluid)  is  very  important  because  only  5/14  "of  the 
work  done  by  buoyancy  appears  as  kinetic  energy  of  mean 
motion,  ...;  the  balance  (9/14)  must  go  into  turbulence 
and  is  eventually  dissipated"  (Ref.  7).  Flow  field  inter¬ 
ference  could  significantly  decrease  the  kinetic  part  of 


93 


the  buoyant  energy  and  allow  the  turbulent  portion  to  domi¬ 
nate;  this  perturbation  in  energy  partitioning  could  cause 
the  strongest  flows  to  take  the  form  of  large  eddies  as 
evidenced  in  the  MBI 1-2  dust  cloud  movies.  This  would 
suggest  that  multiple  burst  cloud  volumes  may  be  determined 
by  turbulent  flow  as  opposed  to  the  strong  vortical  flow 
in  single  burst  dust  clouds. 

The  possible  a, 18  percent  enhancement  in  the  MBI I  dust  cloud 
mass  loading  ratio  may  be  due  entirely  to  the  uniqueness 
of  the  hexagonal  charge  configuration  and  the  simultaneity 
of  the  detonations;  i.e.,  this  particular  configuration 
may  represent  the  "ideal"  one  for  enhanced  flow  fields 
and  mass  lofting  ability.  For  these  reasons,  any  inferred 
mass  loading  enhancement  would  not  have  broad  applicability 
in  "real  life"  attack  scenarios. 
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APPENDIX  A 

FALLOUT  DATA  ANALYSIS 


A-l  INTRODUCTION 

The  procedures  and  preliminary  results  of  the  MBII  fallout 
experiments  fielded  by  SAI  have  been  described  in  reference  1.  Because 
of  anomalously  high  wind,  the  fallout  experiments  were  not  totally 
successful  in  that  fallout  data  was  not  obtained  from  the  single  burst 
cloud,  hence,  early  time  mass  deposition  and  particle  size  distribution 
comparisons  of  single  and  multiple  burst  fallout  data  could  not  be 
accomplished.  Analysis  of  the  fallout  from  the  multiple  burst  cloud 
has  yielded  meaningful  results  concerning  the  source  of  the  early  time 
fallout.  Further,  since  some  fallout  samples  were  obtained  from  regions 
of  the  cloud  where  airborne  sampling  was  accomplished,  a  more  meaning¬ 
ful  qualitative  interpretation  of  the  in-si tu  spectral  data  was 
possible  (since  fallout  data  is  a  time  integrated  sample,  direct  quan¬ 
titative  comparisons  could  not  be  made  with  the  in-si tu  sampling  data). 

A-2  EXPERIMENT  RESULTS 

Fallout  experiments  were  fielded  for  both  events  in  the 
MISERS  BLUFF  II  test  series.  Due  to  unfavorable  and  anomalous  wind 
conditions  existing  in  the  vicinity  of  the  test  bed  for  the  single 
burst  detonation  (23  June  1978),  no  fallout  samples  were  acquired  in 
the  56  fallout  trays  fielded  for  that  particular  test.  During  late 
August  1978,  approximately  70  fallout  trays  were  set  out  for  the  fall¬ 
out  experiment  fielded  for  the  multiple  burst  test.  Figure  A.l  is  a 
topographical  map  of  the  immediate  area  surrounding  ground  zero  and 
depicts  the  locations  of  the  fallout  trays  fielded  for  this  event. 
Favorable  wind  conditions  on  test  day  made  it  possible  to  obtain  18 
dust  particle  samples  from  the  multiple  burst  dust  cloud;  mostly  to  the 
north  and  north-northwest  of  GZ.  The  locations  of  the  samples  are 
shown  in  Figure  A-2.  Each  unit  on  the  grid  equals  500  feet  and  the 
grid  center  (0,  0)  is  the  single  burst  GZ.  The  solid  squares  (  13) 
identify  those  fallout  trays  which  contained  measurable  fallout  while 
the  empty  squares  depict  areas  of  no  recorded  fallout.  The  squares 


Figure  A. 1.  Fallout  Tray  Locations  for  MISERS  BLUFF  II-2  FigureA.2.  Fallout  Sample  Locations  For  The 
Multiple  Burst  Event  MISERS  BLUFF  II-2  Dust  Cloud 


with  diagonals  (6)  mark  the  location  of  dust  samples  that  consisted 
of  small  quantities  of  very  fine  particles  (<  90  um)  probably  lofted 
by  the  blast  wave  passing  over  the  loose  surface  soil.  These  samples 
were  not  attributed  to  dust  cloud  fallout.  The  X's  show  locations  of 
the  fallout  trays  that  did  not  survive  the  high  dynamic  pressures 
within  1000  feet  of  the  "GZ"  for  the  array  of  bursts. 

A-3  FALLOUT  DATA  ANALYSIS 
A-3.1  Sample  Sieving 

In  order  to  facilitate  correlating  the  fallout  data  with  the 
early  time  PSD*  data  obtained  by  aircraft  sampling  the  cloud  stem  and 
cap,  particle  size  distributions  were  obtained  for  the  fallout  samples. 
This  was  accomplished  by  dry  sieving  the  fallout  samples  using  standard 
sieving  meshes  and  plotting  soil  gradation  curves.  These  curves,  which 
yield  a  percentage  of  the  sample  that  is  finer  (or  coarser)  by  weight 
than  the  total  sample  weight,  were  then  used  to  construct  mass  distri¬ 
butions  for  each  of  the  fallout  samples.  The  ten  standard  sieve  open¬ 
ings  used  in  the  fallout  analysis  were  2000,  1400,  1000,  710,  500,  250, 
90,  63,  45,  and  38  microns,  although  not  all  of  the  samples  contained 
particles  whose  sizes  spanned  this  large  range  of  mesh  sizes. 

Of  the  13  dust  cloud  fallout  samples  obtained  from  the  multiple 
burst  test,  only  6  contained  sufficient  dust  (N  0.3g)  to  obtain  the 
accurate  (O.Olg  resolution)  weight  measurements  required  to  generate 
soil  gradation  curves.  The  remaining  7  fallout  samples  were  sieved, 
but  only  to  determine  the  maximum  particle  sizes  (for  both  solid  par¬ 
ticles  and  aggregates)  and  to  qualitatively  estimate  the  range  of 
particle  diameters  that  constituted  the  majority  of  the  particles  in 
the  sample.  Table  A.l  summarizes  the  results  obtained  from  dry  sieving 
the  13  fallout  samples  from  the  MITERS  BLUFF  1 1-2  dust  cloud. 


*Particle  Size  Distribution 
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ons  for  MISERS  BLUFF  1 1 -2  Fallout  Samples 


The  data  contained  in  Table  A.l  for  fallout  samples  1-6  have 
been  plotted  on  standard  soil  gradation  curve  forms  and  are  presented 
in  Figures  A. 3  through  A. 5.  (It  should  be  noted  here  that  our  proce¬ 
dures  used  in  sizing  the  fallout  samples  were  not  the  same  as  those 
used  in  standard  soil  mechanics.  There,  the  aggregates  are  broken 
down  into  their  smallest  component  grain  sizes  prior  to  sieving  and, 
hence,  much  of  the  medium  and  coarse  "sand"  evident  on  the  fallout 
gradation  curves  would  not  be  found  if  standard  procedures  were  used.) 
Visual  analysis  of  the  large  particle  fallout  revealed  that  the  vast 
majority  of  the  particles  n,  250  pm  were  actually  aggregates  that  con¬ 
sisted  of  very  fine  sand  or  silt  cemented  together.  This  is  mentioned 
to  point  out  that  caution  should  be  exercised  if  comparisons  are  made 
between  the  fallout  soil  gradation  curves  and  the  in-si tu  soil  char¬ 
acterization  performed  by  the  U.S.  Army  Waterways  Experiment  Station 
(WES)  for  the  MISERS  BLUFF  II  test  bed  soil. 

A-3.2  Mass  and  Size  Distributions  of  Fallout  Samples 

With  the  soil  gradation  data  obtained  from  dry  sieving  the 
fallout  samples,  differential  curves  approximating  the  mass  distribu¬ 
tions  of  each  sample  were  derived.  Since  the  fallout  experiment  was 
designed  to  obtain  information  on  the  early  time,  large  particle  en¬ 
vironment  within  the  cloud,  emphasis  was  placed  on  particle  sizes 
>  250  pm. 

Considering  the  mesh  sizes  used  in  the  fallout  analysis,  the 
majority  (>_  97  percent  by  weight)  of  the  larger  particles  in  each 
sample  could  be  lumped  by  mass  (AM(i)),  into  six  distinct  size  classes 
or  bins.  The  range  of  particle  diameters  in  each  bin  was  limited  by 
the  mesh  opening  of  the  sieve  in  which  the  particles  were  contained 
(a.j)  and  the  next,  successively  larger  mesh  size  (a^  +  j).  The  med¬ 
ian  particle  size  (a^)  in  each  bin  was  represented  by: 


Figure  A. 3.  MISERS  BLUFF  1 1 -2  Fallout  Sample  Soil  Gradation  Curves 
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Table  A. 2  summarizes  the  range  of  each  size  bin  (i),  the  bin  widths 
(Aa_j),  and  the  median  bin  diameters  (a.)  used  in  the  MISERS  BLUFF  II-2 
fallout  analysis. 


Table  A. 2 

Size  Bins  for  MBII-2  Fallout  Particles 


Bin  (i) 

vai+1  (um) 

Aai  (pm) 

ai  (i 

1 

90-250 

160 

170 

2 

250-500 

250 

375 

3 

500-710 

210 

605 

4 

710-1000 

290 

855 

5 

1000-1400 

400 

1200 

6 

1400-2000 

600 

1700 

Values  of  AM(i)  for  each  size  bin  can  be  obtained  directly  from 
Table  A.l.  These  values  were  obtained  for  the  six  samples  that  con¬ 
tained  a  sufficient  amount  of  dust  for  accurate  weight  determinations 
(Samples  No.  1-6).  Finally,  with  the  values  of  AM(i)  known  for  each 
bin  (i),  differential  mass  distributions  were  generated  for  each  fall¬ 
out  sample  using  the  functional  form  of: 


•  M(  i ) 
'.a . 


which  is  an  approximation  to  the  continuous  function 


From  the  relationship: 


a 


dM 

da 


m(a) 


d_N 

da 


3 

where  m(a)  =  Ka 

and  dN  is  the  number  of  particles  in  da,  we  can  obtain 


d_N 

da 


(1) 


(2) 
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This  continuous  function  can  then  be  approximated  by  the  expression: 

where  . 

rrp 

From  Equation  (3)  it  can  be  seen  that  curves  proportional  to  the  fall¬ 
out  sample  size  distributions  (dN/da)  can  be  approximated  by  plotting 
values  of: 


Figures  A. 6  through  A. 8  contain  plots  of  both  AM(i)/Aa^  and 

AN(i)/Aa..  vs.  a..  for  fallout  sample  numbers  1-6;  the  smooth  curves 

drawn  through  these  points  represent  dM/da  and  dN/da  vs.  a.  Again, 

a  certain  amount  of  caution  should  be  exercised  in  interpreting  the 

- 1  8 

size  distributions  depicted  in  these  figures.  For  example,  K  -  10 
was  chosen  only  to  allow  both  distributions  to  be  plotted  on  the  same 
scale,  therefore,  the  number  of  particles  obtained  by  integrating  the 
size  distribution  between  two  particle  diameters  would  be  a  gross  over¬ 
estimate.  The  curves  do,  however,  approximate  the  shape  (slope)  of 
the  particle  size  distributions  over  the  size  range  of  interest  as¬ 
suming  the  particle  density  (p)  is  not  a  function  of  particle  size. 

In  interpreting  the  mass  distributions,  it  should  also  be 
remembered  that  the  plotted  maximum  values  in  the  distributions  cor¬ 
respond  to  a  range  of  particle  diameters  that  is  represented  by  a 
median  value  (see  Table  A.I.).  Fortunately,  the  mesh  sizes  selected 
for  dry  sieving  the  samples  limit  the  actual  particle  diameter  corres¬ 
ponding  to  the  maximum  value  to  be  off,  at  most,  145  urn  from  the  plotted 
median  diameter. 


AN(i)  _  -1 

A  a  .  "  ^ 
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A-3.3  Source  of  MISERS  BLUFF  1 1 -2  Close-In  Fallout 

In  order  to  determine  the  source  of  the  close-in  fallout  from 
the  MISERS  BLUFF  1 1-2  dust  cloud,  it  was  necessary  to  develop  a  fallout 
model  that  would  adequately  describe  the  gross  features  of  the  f'.I lout 
deposition.  By  calculating  dust  particle  trajectories  based  on  local 
wind  conditions  and  particle  fall  velocities,  and  then  tracinq  these 
trajectories  backwards  in  time  to  possible  source  regions  within  the 
cloud,  it  soon  became  apparent  that  the  cloud  stem  was  the  only  viable 
source  for  the  close-in  fallout. 

A-3.3. 1  Fallout  Modeling 

The  first  step  in  attempting  to  model  the  close-in  fallout  was 
to  see  how,  as  a  function  of  time,  the  cloud  stem  and  cap  traversed 
the  array  of  fallout  trays.  The  dust  cloud  cap  and  stem  track  data 
used  in  the  analysis  are  depicted  in  Figure  A. 9.  These  data  were  ob¬ 
tained  by  Technology  International  Corporation  (TIC)  from  photographic 
analysis  of  the  MISERS  BLUFF  1 1 -2  dust  cloud.  An  obvious,  and  impor¬ 
tant  feature  to  note  is  that  the  cloud  cap  moves  to  the  north  much 
faster  than  the  cloud  stem  from  T=0  to  T+5  minutes;  the  reason  for  this 
being  that,  due  to  height  differences,  the  horizontal  motion  of  each 
portion  of  the  cloud  is  being  controlled  by  '-ifferent  wind  conditions. 

The  cloud  stem  track  data  also  revealed  that  the  low  level 
winds  (surface  to  approximately  2000  feet  AGL)  measured  at  the  Planet 
Ranch  headquarters  at  detonation  time  could  not  be  used  in  a  definitive 
fallout  analysis.  The  low  level,  orographic  effects  of  the  terrain 
in  the  vicinity  of  GZ  (see  Figure  A. 1)  made  it  necessary  to  slightly 
revise  the  measured  winds  to  account  for  the  early  time  stem  motion. 

The  wind  profile  used  in  the  fallout  model  is  given  in  Table  A. 3  and 
was  derived  from  photographed  stem  motion  and  upper  level  winds  mea¬ 
sured  at  detonation  time. 
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CLOUD  STEM 


Figure  A. 9.  MISERS  BLUFF  1 1-2  Dust  Cloud  Stem  (a)  and  Cap  (b)  Movement  From  T+l  to  T+5  Minutes 
After  Detonation 


Table  a.t. 


MISERS  BLUFF  1 1-2  200  Foot  Layer  Averaged  Winds 


A1  titude 

Winds 

A1 ti tude 

Winds 

(KFT  AGL) 

(Deq/Knots) 

(KFT  AGL) 

sfc  - 

IB 

1 

200  - 

1  1  ■ 

2000  - 

2200b 

190/12 

■ 

1  ■  ■ 

2200  - 

2400b 

186/13 

• 

2400  - 

2600b 

182/14 

I 

120/05 

2600  - 

2800b 

178/14 

120/05 

2800  - 

3000b 

713/15 

1 

3000  - 

3200b 

162/12 

aderived  from  stem  track 
b measured  at  time  of  detonation 

In  addition  to  the  wind  data  available,  particle  fall  veloci¬ 
ties  had  to  be  estimated  and  used  in  conjunction  with  wind  data  in 
order  to  calculate  fallout  particle  trajectories.  Unfortunately,  par¬ 
ticle  fall  velocities  will  vary  not  only  with  particle  area  and  mass, 
but  also  from  region  to  region  within  the  cloud,  depending  on  the  flow 
field  being  experienced.  To  keep  the  fallout  model  as  simple  as  pos¬ 
sible,  it  was  assumed  that,  at  any  given  time,  some  particles  of  all 
sizes  within  the  dust  cloud  will  be  falling  at  their  terminal  fall 
velocities  and  that  their  trajectories  vi 1 1  be  determined  by  the 
velocities  and  the  local  wind  conditions.  With  the  trajectories  calcu¬ 
lated  for  various  sized  particles  (see  Appendix  D  for  terminal  fall 
velocity  calculations)  it  was  possible  to  begin  at  a  fallout  tray  loca¬ 
tion  and  trace  the  trajectory  backward  in  time  to  source  regions  within 
the  dust  cloud. 
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Figure  A . 10(a)  shows  the  fallout  sample  locations  and  depicts 
the  trajectories  that  various  sized  particles  would  have  to  follow 
from  the  dust  cloud  to  be  collected  in  several  of  the  fallout  trays. 

The  numbers  at  intervals  along  the  trajectories  indicate  the  altitude 
of  the  particle  at  that  (x,  y)  location  on  the  grid.  The  envelope  of 
all  particle  trajectories  terminating  at  the  sample  locations  should 
define  the  region  of  the  dust  cloud  that  was  the  source  of  the  fallout. 
Figure  A . 10(b)  shows  this  envelope  and  also  the  outlines  of  the  cloud 
stem  and  cloud  cap  widths  from  T+l  to  T+6  minutes.  It  is  quite  evi¬ 
dent  from  this  figure  that  the  cloud  cap,  at  least  that  portion  of  the 
cap  lying  outside  the  outline  of  the  stem,  contributed  virtually  nothing 
to  the  close-in  fallout  deposition. 

To  test  the  hypothesis  that  the  HISERS  BLUFF  II-2  cloud  stem 
was  the  source  of  the  large  particle  fallout  from  the  dust  cloud,  the 
simple  fallout  model  was  applied  to  dust  particles  within  the  stem 
region  of  the  cloud.  Stem  width  and  height  data,  as  a  function  of 
time-after-burst  (TAB),  was  provided  by  TIC  analysis  of  photographic 
data. 

In  the  fallout  model,  the  cloud  stem  at  any  time  is  represented 
by  a  vertically  oriented  cylinder  of  dust,  hence,  projecting  a  circu¬ 
lar  area  in  the  x-y  plane.  In  actuality,  at  times  after  T+2  minutes, 
the  cloud  stem  began  to  slant  as  the  top  of  the  stem  experienced  higher 
wind  velocities,  and  its  projection  in  the  x-y  plane  was  an  ellipse. 

The  elliptical  shape  of  the  stem  projection  would  tend  to  lengthen, 
but  not  broaden,  the  shape  of  the  resulting  fallout  pattern. 

The  vertical  cylinder  of  dust  was  divided  into  200  foot  layers 
so  that  dust  particles  in  each  layer  could  be  subjected  to  the  wind 
profile  existing  at  detonation  time  (TableA.3).  Again  it  was  assumed 
that  at  least  some  of  the  particles  in  each  layer  would  be  falling  at 
their  terminal  fall  velocities.  Dust  particle  trajectories  were  then 
calculated  for  300,  500,  and  1000  urn  particles  based  on  upper  level 
winds  and  their  respective  fall  velocities  (see  Figure  A  . 10(a)) .  Based 
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Fallout  Particles  >  300  ; 


on  the  calculated  trajectories,  particles  of  varying  sizes  in  each 
dust  layer  were  transported  to  the  ground  and  contours  (envelopes) 
were  drawn  on  the  surface  depicting  the  area  in  which  the  single  size 
particles  from  each  layer  would  land. 

Figures  A.ll(a)-(d)  reveal  the  expected  fallout  deposition 
patterns  from  the  cloud  stem  assuming  particle  fall  begain  at  2,  3,  4, 
and  5  minutes  after  burst.  Rg  and  Hs  refer  to  the  stem  radii  and 
heights,  respectively  at  each  time-after-burst.  The  outline  of  the 
cloud  cap  is  also  shown  on  each  one  of  the  figures  to  show  that  early 
time  fallout  from  this  region  of  the  cloud  was  not  detected  in  the 
fallout  trays.  Except  for  predicting  fallout  in  three  trays  where 
no  fallout  was  collected,  and  no  dust  in  one  tray  where  fallout  was 
sampled,  it  appears  that  the  simple  fallout  model,  which  uses  the  cloud 
stem  as  the  only  source  of  fallout  particles,  adequately  describes 
the  gross  features  of  the  MISERS  BLUFF  1 1 - 2  dust  cloud  fallout. 
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APPENDIX  B 


DETAILED  DUST  PARTICLE  SPECTRAL  DATA  DERIVED  FROM 
IN  SITU  SAMPLING  OF  THE  MISERS  BLUFF  II  DUST  CLOUDS 

B-l.  INTRODUCTION 

The  basic  format  of  the  dust  particle  spectral  data  requested 
from  Particle  Measuring  Systems,  Inc.  (PMS)  was  a  two-second  resolution 
time  history  of  particle  number  concentrations  (#/cc),  by  particle  size, 
over  the  size  range  of  particles  anticipated  for  the  MISERS  BLUFF  II 
dust  clouds. 

For  three  of  the  instruments  employed  in  the  sampling  operation 
(FSSP,  LAS-X,  and  0AP-200X),  obtaining  the  two-second  (%400  ft,  -120m) 
resolution  within  the  dust  clouds  was  possible  via  the  instrument's 
interfacing  with  the  on-board  Data  Acquisition  System  (DAS).  (For  a 
more  detailed  discussion  on  the  instrumentation  employed  in  the  in-situ 
sampling  operation,  refer  to  Reference  2  ) .  Figure  B.l  is  an  example  of 
the  spectral  data  derived  from  the  FSSP  probe  output.  Bin  numbers 
106-120  correlate  with  median  particle  diameters  of  3-45  pm  in  3  pm 
increments,  i.e.,  1.5-4. 5  pm,  4. 5-7.5,  ...,  43.5-46.5  pm.  Time  of 
day  (in  seconds)  is  given  in  ten-second  increments  along  the  bottom  of 
the  figure;  cloud  entry  and  exit  times  are  depicted  by  arrows.  Each 
vertical  column  of  numbersrepresents  the  particle  number  concentration 
(log  (f/cm  ))in  each  size  class  of  particles  every  two  seconds.  Two- 
second  totals  (log  values)  appear  just  above  the  time  of  day. 

Figure  B.2  is  similar  spectral  data  derived  from  the  Two 
Dimensional  Optical  Array  Imaging  Probe  (0AP-2D-C).  Bin  numbers  1-20 
correlate  with  median  particle  diameters  of  30-500  ,im  in  25  ;.m 
increments.  Spectral  data  from  the  0AP-2D  probes  are  recorded  only 
when  particles  are  present  and,  unlike  the  three  probes  mentioned  above, 
this  digital  image  data  is  stored  in  a  buffer  prior  to  being  recorded. 
Since  a  buffer  is  "dumped"  on  tape  only  when  it  becomes  full,  it  is 
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Pass  7  Data  From  First  Pass 
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possible  to  acquire  spectral  data  at  the  end  of  one  sampling  pass  and 
have  this  data  dumped  during  the  next  pass  when  the  buffer  fills.  This 
delay  is  also  true  within  regions  of  the  cloud  where  large  particle 
concentrations  are  low. 

For  this  reason,  obtaining  the  required  two-second  data 
resolution  within  the  MB  1 1  dust  clouds  required  a  significant  software 
development  effort  by  RMS.  Along  with  each  image,  the  elapsed  time 
since  the  last  particle  encounter  is  also  stored  in  the  buffer,  and  these 
times,  along  with  the  buffer  dump  times,  are  recorded.  Hence,  it  was 
possible  to  "back  out"  the  "location"  (time  of  acquisition)  of  each 
particle  encounter  and  place  this  data  in  its  correct  two-second 
interval.  One  remaining  problem  was  that  if  more  than  approximately 
7  seconds  (0AP-2D-C)  elapsed  between  particle  encounters,  the  instru¬ 
ment's  internal  "clock"  may  have  recycled  several  times  and  the 
"location"  of  these  particles  would  be  ambiguous  by  some  multiple  of 
7  seconds.  Fortunately,  by  comparing  these  "questionable"  data 
points  with  the  FSSP  spectral  data,  it  was  possible  to  determine 
(judgemental ly)  the  most  probable  time  of  2D  data  acquisition.  (The 
0AP-2D-P  instrument's  internal  "clock"  recycled  every  ^28  seconds). 

Figure  B.2  is  a  graphic  example  of  how  data  was  acquired  during 
Pass  6  in  the  MBII-1  dust  cloud  and  was  not  recorded  until  the 
buffer  was  filled  and  dumped  during  Pass  7.  (The  actual  time 
between  images  in  this  example  is  98  seconds,  i.e.,  14  clock  cycles). 

This  effort  was  successfully  completed  and  we  feel  confident  that  the 
MISERS  BLUFF  II  dust  cloud  spectral  data  represents  the  most  detailed 
set  of  data  made  to  date  within  a_ny^  cloud  -  dust  or  water  -  with  the 
2D  instruments. 

B-2 .  DETAILED  DUST  PARTICLE  SPECTRAL  DATA  PRESENTATION 

Tables  B.l  and  B.2  are  computer  listings  of  all  the  MISERS 
BLUFF  II  two-second  spectral  data  that  were  used  in  the  analyses 
contained  in  this  report.  Table  B.l  contains  the  spectral  data 
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acquired  during  sampling  passes  1-15  in  the  MISERS  BLUFF  II  single 
burst  dust  cloud  (MBII-1)  while  Table  B.2  contains  similar  data  from 
passes  1-11  in  the  multiple  burst  dust  cloud  (MBII-2).  The  format  for 
both  tables  is  basically  the  same;  the  only  difference  being  Table  B.2 
contains  spectral  data  for  two  more  particle  size  ranges  than  does 
Table  B.l.  This  is  due  to  the  0AP-2D-P  instrument  malfunctioning 
during  the  MBII-1  sampling  mission  and  the  addition  of  the  LAS-X  to 
the  instrumentation  suite  fielded  for  the  MBII-2  sampling  mission 
(see  Reference  1). 

In  both  tables,  the  time  in  seconds  after  detonation,  is  given 
in  the  extreme  left  and  right  hand  columns  of  each  page.  The  asterisks 
refer  to  data  intervals  where  the  2D  spectral  data  has  been  added  to 
a  sampling  pass  as  explained  in  section  B-l.  Extinction  coefficients 
and  radar  reflectivity  factors  corresponding  to  the  asterisked  data 
intervals  have  not  been  corrected  and  are  lined  out. 

Four  particle  parameters:  (1)  particle  number  concentration, 

(2)  mass  concentration,  (3)  extinction  coefficient,  and  (4)  radar 
reflectivity  factor  are  listed  for  each  two-second  data  acquisition 
interval.  Of  the  four  parameters,  only  the  particle  number  concen¬ 
tration  was  measured;  all  other  parameters  have  been  calculated  from 
these  measured  values  assuming  spherical  particles. 

B-2.1  Calculation  of  Mass  Concentrations  (g/m  ) 

As  previously  mentioned,  the  only  measured  quantity  available 
from  the  PMS  data  tapes  was  the  particle  number  concentration  (PNC) 
for  each  size  class  of  particles  during  each  two  second  data  acquisition 
period.  Tor  each  size  class,  the  mass  concentration  (MC)  is  given  by 

MC(n/m3 )  -  PNC ( « / cm3 )  1 06 ( cm3 /in3 )  Particle  Mass  (g) 
or 

3 

MC(g/m3)  -  (PNC) (  ).  x  106 

where  (I  is  tin'  particle  diameter  in  cm  and 

.  3 

in  q/cm  . 
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is  the  particle  density 


For  each  two-second  interval,  these  values  are  summed  over 
all  size  classes  (44  for  MBII-1  data  and  64  for  MBII-2  data)  to  yield 
the  total  mass  concentration,  whose  logrithmic  value  is  tabulated  in 
Tables  B.l  and  B.2.  The  percent  contribution  to  the  total  mass  con¬ 
centration  from  several  particle  size  ranges  has  also  been  calculated 
for  tabulation. 

B-2.2  Calculation  of  Extinction  Coefficients  (1/km) 

The  extinction  coefficient  (EC)  for  each  particle  size  class 
in  a  two-second  data  interval  has  been  approximated  by  the  Hie 
formula : 

EC(l/km)  =  2[PNC(#/cm3)103(cm/km)  Particle  Area  (cm3)] 
or 

...  C 

EC(l/km)  =  (PNC) (4?-)  x  1(T 

where  d  is  the  particle  diameter  in  cm  (this  significantly  over¬ 
estimates  for  5d  .  optical  wavelength). 

These  values  are  also  summed  over  all  particle  size  classes 
present  in  each  two-second  data  interval  to  obtain  the  total  value 
whose  logrithmic  values  have  been  tabulated.  As  with  the  mass 
concentrations,  the  percent  contribution  from  several  particle  size 
ranges  has  also  been  tabulated. 

6  3 

B-2.3  Calculation  of  Radar  Reflectivity  Factors  (mnr/nr) 

Radar  Reflectivity  Factors  (RRF)  have  also  been  calculated 
and  tabulated  for  each  two-second  data  acquisition  interval.  For 
each  particle  size  class  in  a  data  acquisition  interval,  the  RRF 
has  been  approximated  by  the  Rayleigh  formula: 

RRF(mm6/m3)  =  PNC(«/cm3)106(cm3/m3)  10('(mm6/cm6 )  (Part  i cle  Diameter 
or 

RRF (mm6 /m3)  -  PNC(d6)xl013 
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where  d  is  again  the  particle  diameter  in  cm  (this  significantly  over¬ 
estimates  for  5d  >  radar  wavelength). 

The  total  radar  reflectivity  factor  in  each  two-second  data 
interval  is  obtained  by  summing  over  all  particle  size  classes  present 
in  that  interval.  Logrithmic  values  of  total  radar  reflectivity  factor, 
plus  the  contribution  to  this  total  by  several  particle  size  ranges, 
are  tabulated  in  Tables  B.l  and  B.2. 
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APPENDIX  C 


THE  MISERS  BLUFF  1 1-2  AEROSOL  ENVIRONMENT 
C-l  INTRODUCTION 

The  presentation  of  spectral  data  obtained  from  the  LAS-X 
instrument  (particle  size  range  0.1  to  6  pm)  is  being  included  as  an 
appendix  to  this  document  since  these  particle  measurements  were  not 
common  to  both  the  single  and  multiple  burst  dust  clouds.  The 
addition  of  this  particular  probe  to  the  in-situ  sampling  instrumentation 
suite  was  prompted  by  an  unexpected  wavelength  sensitivity  in  lidar 
backscattering  measurements  obtained  by  SRI  International  during  the 
single  burst  test.  It  is  hoped  that  the  inclusion  of  this  spectral 
data  obtained  from  the  MBII-2  dust  cloud  will  aid  in  interpreting  the 
lidar  backscatter  measurements  made  at  10.6  pm,  1.06  ;.m  and  0.532  pm 
during  both  events.  This  may  be  possible  since  the  very  small  particle 
environment  was  probably  similar  in  both  dust  cloud  environments  due 
to  the  ANFO  and  silt  commonality. 

C-2  PARTICLE  SIZE  DISTRIBUTIONS 

To  be  consistent  with  Sections  2  and  3  of  this  report,  t he 
size  distributions  obtained  for  dust  particles  _6  ,.m  will  be  presented 
only  for  those  sampling  passes  that  were  made  while  ascending  to  the 
cloud  top,  i.e.,  passes  2-11  in  the  MBII-2  dust  cloud.  Further,  Lite 
LAS-X  output  was  not  recorded  during  sampling  passes  4,  5,  and  11  so 
pass-averaged  size  distributions  are  not.  available  for  those  three  passes. 

Figures  C-l  through  C-7  depict  the  average  size  distribution 
of  particles  measured  by  the  LAS-X  instrument  during  ascending, 

MBII-2  sampling  passes.  The  solid  circles  depict  the  particle  size 
distributions  of  particles  5-25  ,.m  in  diameter,  as  measured  by  the 
I SSP  instrument,  and  are  included  only  to  show  the  excellent  agreement 
between  the  two  probes.  (For  various  reasons,  the  solid  circles  at, 

5  ,,m  are  less  reliable  -  see  Reference  2). 
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Fiqure  C.j.  MISERS  BLUET  1 1 -2 ,  Pass  8,  Ae-osol  Particle 
Size  Distribution 
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Fiqure  C.7.  MISFRS  BL UFT  11-2 ,  Pass  ID,  Aerosol  Particle 
Si zp  Distribution 


For  anyone  who  is  focusing  on  dust  particle  sizes  3  ,.m  in 
their  analysis  of  the  MISERS  BLUFF  1 1-2  dust  cloud,  particle  number 
concentrations  (f/cc)  and  calculated  mass  concentrations  (g/m  ),  in 
two-second  data  intervals  ('^400  foot  resolution)  within  the  cloud, 
can  be  easily  extracted  from  Appendix  B  (Table  B.2)  to  supplement  the 
size  distribution  data  presented  in  this  appendix. 

Since  the  acquisition  of  spectral  data  for  particle  sizes  less 
than  3  pm  was  not  required  for  cloud  mass  loading  comparisons,  only 
general  comments  will  be  made  concerning  these  size  distributions. 

First,  it  appears  that  between  0.1  and  6  urn,  the  size  distribution 
of  these  particles  cannot  be  adequately  described  as  either  exponential 
or  power  law  in  nature.  This  fact  could  affect  an  analysis  of  lidar 
measurements  where  an  assumed  particle  size  distribution  was  used  to 
calculate  backscatter  coefficients.  Second,  there  appears  to  be  a 
"flattening"  of  the  size  distributions  between  "•0.2  pm  and  • 0.7  pm; 
in  fact,  most  of  the  pass-a veraged  distributions  indicate  a  "hump" 
in  the  distribution  at  -0.3  pm.  This  could  lead  to  a  significant 
perturbation  in  any  analysis  that  assumed  the  particle  size  distribution 
could  be  represented  by  a  monotonical ly  decreasing  function. 


APPENDIX  D 

DUST  PARTICLE  TERMINAL  FALL  VELOCITIES 


Since  much  of  the  analyses  of  the  MISERS  BLUFF  II  dust  clouds 
relied  heavily  upon  assumptions  concerning  dust  particle  terminal  fall 
velocities,  it  was  felt  a  brief  appendix  pertaining  to  particle  fall 
velocities  was  warranted.  Although  all  analyses,  for  lack  of  more 
definitive  information  on  dust  particle  shapes  over  the  size  range 
of  interest, were  based  on  "rough"  spherical  particles,  this  appendix 
also  contains  terminal  fall  velocities  of  circular  discs,  a  shape  not 
totally  inconsistent  with  the  "mica  flakes"  reportedly  found  by  PMS 
in  thei  r ai rcraft  engine  filter. 

U-l  SPHERICAL  PARTICLES 


For  a  spherical  particle  of  mass  m  falling  at  its  terminal 
velocity  vt  through  air  of  density  pa  ,  the  force  balance  is 

«d3g  MaVt2CDAc 

mg  =  -5—-  = - 2 -  (1) 


where  a  is  the  cross-sectional  area  of  a  particle  with  diameter  d, 

(>  is  the  particle  density,  and  C'D  is  the  particle's  drag  coefficient. 
The  drag  coefficient  is  an  empirical  function  of  the  Reynolds  number 
Re  and  is  given  by 


CD  =  Cj  +  (C2/Re)  (2) 

where  C,  and  C0  have  been  experimental  1 y  determined  for  various 
12 

shaped  particles.  For  "rough"  spheres  Cj  =  0.6  and  C^  =  36.  By 
comparison,  for  smooth  spheres,  C^  =  0.4  and  C^  =  24,  which  indicates 
an  expected  decrease  in  drag  coefficient  for  a  niven  Reynolds  number. 
The  Reynolds  number  of  a  particle  with  diameter  <1  is  defined  as 

Re  =  ,-a  Vt  d  /  (3) 

where  ,,  is  the  dynamic  viscosity  of  air.  Substituting  (2)  and  (3) 


D-l 


into  (1)  and  letting  Ac  =  ndV4,  we  obtain  a  quadratic  in  which  is 
given  by 

O  c?MVt  4p  dg 

r  X,  £  .  t  P  _  n  /  \ 


it  pad  3pa 


Solving  (3a)  for  the  non-negative  root  of  1/  we  have 


-c2m  +  ^ (C2,i)2  +  ]f  (PpPggCjd3) 


2pacid 


Equation  (3a)  can  be  rewritten  as 


2  ^2^t 

r  v  c  +  Jz _ F_ 

Vt  Re 


4«pdg 


which  indicates  that  for  large  values  of  Re,(i.e.»  large  particles), 
the  first  term  in  2  will  dominate  the  equation  and  v'fc  can  be 
approximated  by 

Vt  r  yj  4,  pdg/3Cl(  a  ,  d  >  0. 1  cm 

Similarly,  for  small  values  of  Re,(i.e.,  small  particles) ,  the  second 
vV  term  will  dominate  the  equation  and  Vf  will  be  approximated  by 


Vt  =  4'  pd  a/3C2"  ’  d  ~  0>01  Cl" 

From  this  equation  we  note  that,  for  small  particles,  the  density 
of  air  through  which  they  are  falling  has  no  direct  effect  on  t he 
particles  fall  velocity. 

From  the  meteorological  data  obtained  during  Phase  1  of  the 
MISERS  BLUFF  II  Cloud  Sampling  Program  (see  Ref  1),  the  air  density 
for  both  events  as  a  function  of  altitude  at  the  test  site  is  given 

i 

approximately  by 


pa  (g/cm3)  =  1.128  x  10"3  exp  (-0.02657Z) 

where  Z  is  the  altitude  in  kft  MSL.  The  dynamic  viscosity  of  air  is 
given  approximately  by 

Pa  (g/cnTsec)  =  (178. 7-Z)xlO~6 

where,  again,  z  is  the  altitude  in  kft  MSL. 

For  the  terminal  fall  velocity  calculations  performed  in  con¬ 
nection  with  the  T+10  minute  and  T+20  minute  cloud  reconstructions, 
an  average  altitude  of  5  kft  was  chosen  as  being  a  "representative" 
altitude  at  which  to  compute  particle  fall  velocities.  The  maximum 
errors  induced  by  this  assumption  over  the  altitude  range  of  interest 
are  typically  <  10  '  for  the  larger  particles  (>1  mm)  and  much  less 
for  the  smaller  particles  (<500  pm).  For  the  calculations  used  in 
the  fallout  data  analysis,  the  "representative"  altitude  chosen  was 
2  kft  which,  over  the  range  of  cloud  stem  heights,  induced  much  smaller 
errors . 

Another  constant  used  in  the  terminal  fall  velocity  calcula¬ 
tions  was: 

g  =  981  cm/sec 

Substituting  the  constants  into  (4)  we  have,  for  "rough"  spherical 
particles , 

-36;:  +  V1250  +  3136  ■  n'  Ad  ^  (5) 

V t  (cm/sec)  = - y/Td . ^ - 

d 

where  0  is  in  cm,.  All  calculations  in  this  renort  assumed  a  value 
of  2.3  g/cm3  for  .  . 

p 

0-2  CIRCULAR  DISCS 

Much  of  the  information  presented  concerning  the  fall  velocities 
of  circular  discs  has  been  extracted  from  Reference  3.  Solving  (1) 
for  Cp  and  multiplying  by  the  square  of  (3)  we  obtain 


X 


(6) 


CD(Re)2 


2md2gp 


V 


x  is  defined  as  the  Best  number  (after  a  person  who  did  much  work  in 
this  area  of  physics)  and  note  that  it  depends  only  on  the  particle 
and  its  position  in  the  atmosphere  -  not  on  speed  and  drag.  For  a 
circular  disc  falling  in  a  horizontal  attitude, 

m  =  t'Q  n  d2h/4  and  Ac  =  n  d2/4  (7) 


where  h  is  the  thickness  of  the  disc,  and  hence. 


X 


2i>  hd2i!  g  /  ti2 

P  3 


(3) 


Using  the  same  constants  as  used  for  the  spherical  particles,  and 
defining  the  aspect  ratio  r  of  a  circular  disc  as  being  the  ratio 
of  its  diameter  to  its  thickness,  we  obtain 


X  =  (4 . 5 lx 103  ,  /  u2)(d3/r) 

a 


(9) 


For  particles  of  various  shapes,  X  and  c  have  been  experimentally 
determined,  and  hence,  from  (6),  the  Reynolds  number  can  be  determined 
as  a  function  of  the  Best  number.  Over  the  limited  range  of  Best 
numbers,  20  1  X  <  2000  and  drag  coefficients  2  <  Cp  <  200,  the 
Reynolds  number  for  a  circular  disc  can  be  approximated  by 

Re  =  0.1  0-786  (10) 


Substituting  (9)  into  (10)  we  obtain  from  (3) 


V 


t 


(cm/sec)  =  74.5 


-.0 . 786 


(ID 


which  approximates  the  terminal  fall  velocity  of  a  circulai  disc. 


i)-4 


Terminal  fall  velocities  for  "rough"  spheres  and  circular 
discs  with  several  aspect  ratios  have  been  calculated  at  an  altitude 
of  5  kft  MSL  and  are  presented  in  Figure  D-l.  As  would  be  expected, 
the  circular  discs  fall  significantly  slower  than  do  the  "rough" 
spheroids  with  the  same  diameter  and  this  difference  increases  with 
increasing  aspect  ratios  of  the  circular  discs.  (In  actuality,  the 
curves  depicting  the  fall  velocities  for  circular  discs  are  slightly 
curved,  but  due  to  the  Reynolds  number  being  approximated  by  (10), 
they  appear  as  straight  lines  over  the  range  of  particle  diameters 
that  the  approximation  is  valid.) 
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